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= is equal by definition to 

sgn(x) sign of x 

A\B elements of set A not in B 

A ± positive or negative numbers 

of A 

A x H cartesian product of two sets 

$t(z), ^s(z) real, imaginary part of z G C 
A = B A and IB are two equivalent 

matrices 



0{ Xl ,...,x n }(k) (set of) smooth function(s) of 
x = {xi, . . . , x n } whose deriva- 
tives up to order k — 1 vanish 
at {0, ...,0}. 
o, a, p, etc. fixed points or periodic orbits 
A,B,R,A 2 x 2-matrices 
sp(U) spectrum of matrix U 

e* \e\ if I = 3, tJ\e\ otherwise 
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1 Introduction 



Semiclassical methods have had many applications in the description of the quantum mechanical 
behavior of particles whose classical analog is chaotic. These methods are essentially based on the 
knowledge of the classical periodic orbits. Aside from these classical contributions, there are other 
ones, usually associated to non-classical effects such as diffraction or tunneling. The main motivation 
of this work is to better understand a particular type of contribution to trace formulae, i.e. the 
complex periodic orbits. We are thus interested in semiclassical expansions which generically take the 
form of an asymptotic series 

Q (h) ~ j2 A iP^ wm/H (i +£°*M R *) • (1.0.1) 

Here, the first sum concerns (periodic) orbits and fixed points of the classical system. Both A[p] 
and W[p] are complex functions of geometrical quantities attached to each p. The prefactor A[p] 
is well defined if the periodic orbits are assumed to be isolated. This is typically the case for fully 
chaotic systems, as for example the Sinai billiard or the Bunimovich stadium. However, no smooth 
potential is known to be fully chaotic and the generic scenario is a coexistence of regular and irregular 
motion in phase space (the so called mixed dynamics). In this case the phase space structure is much 
more complicated than in the fully chaotic regime because the frontier between regular and irregular 
motion is not smooth, but fractal, and its details are very sensitive to changes of external parameters 
controlling the dynamics. 

At the level of periodic orbits, this complexity is reflected on the one hand in a rich structure of 
coexisting stable and unstable orbits and on the other hand by the occurrence of many bifurcations as 
some parameter is varied. Close to a bifurcation several periodic orbits can lie very close to each other, 
and exactly at the bifurcation they coalesce. Semiclassically, the distance between periodic orbits is 
measured in terms of a phase space length-scale fixed by the value of h. For distances smaller than this 



scale, the prefactor A[p] in (1.0.1) diverges, and some uniform approximation treating simultaneously 
the contribution of several orbits has to be implemented. For two-freedom systems, for which a 
complete classification of generic bifurcations exists p4| (see also table I below), this was done in 
Refs. |Iq| , |4q| . Because most bifurcations imply the creation of new or the disappearance of existent 
periodic orbits as a parameter is varied, uniform approximations are not only a good way to solve the 
problem of the divergence of the prefactor, they also take care of the (unphysical) discontinuities that 
would have occurred in fll.O.lj ) if suddenly a term is added (or subtracted) to the interferent sum. In 



fact, the most natural way to obtain a continuous physical picture through the bifurcation is to extend 
the classical dynamics (i.e. the coordinates and momenta) into the complex plane. Then it is easy to 
see, as we will show below, that classically there is no discontinuity (i.e., no creation or destruction) 
at the bifurcation but rather an exchange of periodic orbits flowing from the complex to the real plane 
and vice versa as a parameter is varied. (Complex periodic orbits and more generally complex classical 
trajectories are solutions of the classical equations of motion whose coordinate and/or momenta are 
complex variables. Time is always considered as a real variable). Exactly at the bifurcation point 
several orbits coming from the complex and/or the real plane coalesce. 

The inclusion of complex orbits in semiclassical formulae is thus an important issue and may be 
considered as part of a general program on exact semiclassical mechanics. The aim of this program 
is to compute exact quantum results from semiclassical analysis (see for example the review by A. 
Voros [^] . As Balian and Bloch have pointed out , quantum mechanics can be fully reconstructed 
from classical trajectories if real as well as complex trajectories are included. This program, which for 
generic mixed systems remains for the moment at a formal level, has however been pursued successfully 
for one-dimensional autonomous systems |5^, |54| as for the free motion of a particle on a compact 
surface of constant negative curvature [15| . In spite of the fact that the latter is a prototype of fully 
chaotic motion, the accomplishment of the program in this case heavily relies in the fact that the 
semiclassical expressions are exact. 

The appearance of complex trajectories in exact semiclassical expressions is of course easy to un- 
derstand. Consider for example the tunneling effect. By definition this effect is classically forbidden 
and hence is not taken into account in expressions such as Eq. ( 1.0. 1| ) where only real classical tra- 
jectories are included. It is however remarkable that only complex classical paths (and not arbitrary 



complex paths) are needed in order to correctly describe it ||]. Ref. [31|, where complex periodic 
orbits are used in order to include tunneling effects in the computation of the density of states of a 
mixed system, provides an illustration of this point. Other interesting recent works exploring different 
aspects of tunneling in the mixed or fully chaotic regimes can be found in Refs.||, ||, g ||, |l|, ||. 

As mentioned before, the amplitude ^4[p] in Eq. (1.0.1) depends on the geometry of the classical 
flow and is connected to the linearized motion around p. If p is isolated, it can be expressed in terms 
of the eigenvalues of the monodromy matrix M at p which is, by definition, the linearized part of 
the transversal map <I> after one return in the neighborhood of any point of p. In order to find 
these eigenvalues, one can, for example, diagonalise the quadratic part of the classical Hamiltonian by 



5 



constructing a suitable symplectic chart in which the later is made up of only the relevant geometric 
monomials. This procedure can be generalized to higher orders: it is the object of normal form theory. 
In hamiltonian physics, it has been introduced by El] and ]|] . 

In order to obtain geometrical semiclassical h expansions, one must therefore classify the generic 
situations which may happen in the neighborhood of periodic orbits as an external parameter is 



varied. In [34] |36|] such a classification was performed for one-parameter hamiltonian systems with 
two degrees of freedom. In the original paper, the aim was to prove the genericity in the mathematical 
sense. In a more recent work f37|, the authors use the Lie transformation theory developed in the 



context of many body classical problems [27, ^Cj, and references therein], which is a powerful method 
for obtaining normal forms. It has also the great advantage of giving a systematic method that can be 
numerically implemented. However, because our final aim are semiclassical expansions, the classical 
objects we are interested in are the generating functions - the phases in Feynman integrals - rather 
than the explicit maps themselves. Especially when the bifurcation is of an extremal or transitional 
type (see section |] below), the explicit expression for the generating functions is not trivially obtained 
from Meyer's work. We will give a precise expression for them in section ^. 

Before, the aim of section ^ is to introduce the one-parameter hamiltonian systems with two degrees 
of freedom, and of sections |3| and |] is to give an original presentation of Meyer's classification. We 
will neither use the language of singularity theory nor that of Lie transformations; we will only need 
some elementary algebraic operations and refer to a physicists' intuitive approach of genericity. We 
will systematically use Taylor expansions and make generic statements about the coefficients. 

In section ||, we make a detailed study of the structure of the geometrical classical skeleton based 
on the previous results and include systematically the complex fixed points in the analysis of each type 
of bifurcation. The reader familiar with the theory of normal forms or not interested in the proofs can 
directly go to this section were the essential ingredients are summarized. 

The quantum expansions based on these classical structures are finally introduced in section 
We go through all the different types of bifurcations encountered in two-dimensional systems and 
systematically evaluate the contribution of the real as well as the complex periodic orbits involved. 
The main difficulty is associated with the complex orbits. When one uses the steepest descent method 
in semiclassical approximations one must take into account global properties of the integration domain 
to find the possible deformation of the integration path that determines the contribution of the orbit. 
In particular, the direction in which one crosses the saddle points must be known in order to correctly 
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compute the Maslov indices of the orbits. When the dynamics is embedded in a smooth family, 
the choice of critical points may change discontinuously as the external parameters are varied. This 
produces the well known Stokes'phenomenon extensively studied by Berry & Howls || (see 
also |2q| ). The difficulties when dealing with this highly non-local properties may be attenuated if 
one considers only complex orbits near a bifurcation. Aside from these considerations, our interest 
is to compute the relative weight of the contribution of complex periodic orbits. Several previous 



studies in the mixed regime [30, 31] have shown that some complex orbits give an exponentially small 
contribution. Nevertheless we will argue in section [5] that these exponentially small contributions 
occur in only one type of bifurcation: the extremal case. For all the others the contribution of 
complex periodic orbits - if they exist - is proportional to a power law in h. In other words, these 
complex periodic orbits contribute to akh k -\\ke terms in Eq.( 1.0.l| ). 
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2 Poincare reduction 



Let us consider an autonomous hamiltonian system with two degrees of freedom. Denote G its hamil- 
tonian, denned on a four-dimensional phase space. Let p be a periodic orbit of period T > lying in 
the three-dimensional energy shell H(E). Let us construct another three-dimensional submanifold & 
which is sufficiently small to intersect p only once. Let o = © n p and s(E) = 6f1 E(E) (see figure Of) . 




<f 



Figure 1: Poincare reduction in the neighborhood of a periodic orhit p. 



Let us call $ the angle coordinate defined in a tubular neighborhood of p. On p, it represents 
the arclength coordinate modulo the total length r of p counted in a suitable unit system. Moreover 
there is, by construction, a one-to-one relation between and the time t since everywhere locally 
around p, d'd/dt ^ 0. One can thus consider t to be a smooth function of £(i9) . Let I be the coordi- 
nate canonically conjugate to •& and (p, q) the other two symplectic coordinates. The hamiltonian G 
is locally a function of these variables: G(I,p, q). Since dG/dl = d'd/dt ^ 0, one can make use of 
the implicit function theorem to define a function H (p, q; E) by 

G(I,p, 0,q) = E <=> 1 = Hip, q; tf; E) . (2.0.2) 

From the expression of the partial derivative of the implicit functions we get 
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dp _ dp/dt _ dG/dq dH 

d& ~ ddjdt ~ ~dG/dI ~ + ~dq~ ' ( ' 

dq = dg/dt _ dG/dp _ dH 

d& d&/dt dG/dl dp ' 

Therefore, this construction - due to |55|, chapter XII, §141], see also [R chapter VI, §3] - reduces 
locally the original two freedom system to an equivalent one-dimensional system. The price to pay 
is that the new system is nonautonomous since the new hamiltonian H depends on the new time i?. 
Nevertheless, the reduced dynamics, now living on s(E), is r-periodic due to the periodicity of o is 
an equilibrium point of H obeying, V 

OH dH 

-Q- P =o = and — p=0 = , (2.0.4) 

U r q=0 Uq q=0 

if (p, q) is centered on o . The reduced dynamics mimics the transversal motion stroboscopically with 
period r. The corresponding one-step iteration map &(E; r) defined on s(E) is the so-called Poincare 
map. Exploring what happens around p when we go outside X(-E') is equivalent to considering that the 
reduced dynamics is embedded in a one-parameter {E) unfolding. Since the symplectic chart on s(E) 
is defined in a neighborhood of it, we will adopt a passive point of view by considering that we stay 
for all E in the same bidimensional phase space . 



3 Normal form for an individual hamiltonian 

From the preceding arguments, changing slightly the notation, let us start from a smooth Hamilto- 
nian H(p, q; t) defined on a bidimensional phase space & and which is r-periodic in time (r > 0). We 
will assume that the origin is an isolated fixed point. 

If we denote w = f (g) and V„, = (^ p ) , the equations of motion can be written in the compact 
form 

w(t) = JV w H(w(t);t) (3.0.5) 

where " ' " d = d/dt and J = f ( ° ~\ ) . 
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3.1 Elimination of the time-dependence in the linearized motion 

The classification of the dynamics in the neighborhood of an isolated fixed point is based on the 
linearization of the motion around it. 

Linearizing equation ( |3.0.5| ) around the origin it follows 

w(t) = A(t) w(t) (3.1.1) 

where 

(d p 2 p H(w;t) d 2 q H(w;t)\ 



def 

- 



A(t) = -J Hess [#(«;;*)] 



(3.1.2) 



\d q 2 p H(w;t) d q 2 q H(w;t)J 



-0 



It is natural to expect that the linearized dynamics ( 3.1. 1| ) is roughly akin to the true dynamics 



provided we are close enough to the origin. We will come back to this point at the end of the 
subsection. 

Recall that a symplectic (complex) matrix § is defined by the condition 

t SJJ§ = J (*S denotes the transpose of 8) (3.1.3) 
In two dimensions, this is equivalent to 

detS = l. (3.1.4) 

By definition Q definition 3.1.13] a complex matrix IK is an infinitesimal symplectic matrix if it 
verifies 

JK = -*KJ. (3.1.5) 

Then it can easily be checked that A(i) is an infinitesimal symplectic matrix. As its name suggests, 
the exponential of an infinitesimal symplectic matrix is symplectic. 

Let U(i) be the evolution matrix associated with A(i). In other words, w(t) = V(t)wo is the 
unique solution of equation fl3.1,l|) such that w(0) = wq. We have 

U(t) = A(t)U(t) with U(0) = l= f (J I) . (3.1.6) 

One often writes the formal solution of the latter equation as a time ordered product 

U(t) =Te /o<A(m ' . (3.1.7) 
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It is straightforward to show that U(i) is symplectic for all t. Let L be any time-independent infini- 
tesimal symplectic matrix. Then, being the product of two symplectic matrices, F(t) = e* [U(i)] 1 
is also symplectic. If we make the symplectic, time-dependent change of coordinates w = F(i) w, then 
equation ( |3.1.l| ) becomes w(t) = hw(t). All the information of the original linearized dynamics has 
been hidden in the change of variables. Let us keep in mind that H (p, q; t) is a reduced hamiltonian and 
hence that is a (local) Poincare surface of section centered at one point of an isolated periodic orbit p 
of a higher dimensional system. Then, the time-dependence of the flow defined by equation ( |3.0.5| ) 
can be seen as a transversal projection of the non-reduced autonomous flow as we move along p. The 
time periodicity of the reduced flow comes from the fact that after one cycle along p we are back 
on & '. Hence we should consider only real r-periodic symplectic changes of coordinates, F(r) = i. 
L must then be chosen such that e T ^ = M = U(r), the so-called monodromy matrix. rL exists and 



can be any logarithm of M since M has no zero eigenvalues as a consequence of equation ( |3.1.4 ). 



Nevertheless, it may happen that no real logarithm of M exists. In this case one possibility is to relax 
the constraint on F(i) by retaining only the intersection with the Poincare surface of section every 
second crossing. That means that F(t) must only be 2r-periodic. Because U(2r) = [U(r)] 2 and using 
the fact that every real matrix that has a square root has a real logarithm f37| , theorem 2 of II. E], we 
can indeed choose the infinitesimal symplectic matrix L to be real without loss of generality. This is 
the statement of the Floquet theorem [23] (see also the work of Cherry |L7], [u], [l8| ). There exists 



a symplectic coordinate system in which the linearized motion ( 3.1.1[ ) is governed by a constant real 



infinitesimal symplectic matrix L. Equivalently, in a well chosen symplectic chart, the hamiltonian 
can be written 

H(p,q;t) = ^cp 2 + dpq+-eq 2 + h^(p,q;t) (3.1.8) 



def 



h^(p,q) 



in a neighborhood of the origin after having subtracted the irrelevant value of if(0,0;i). (c, d, e) 6 
[R 3 and is a r-periodic function whose first and second derivatives with respect to p or q vanish 

at the origin for all t. In this case A(t) is a constant matrix A and, from definition ( gXp 




L = A= , (3.1.9) 
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and the monodromy matrix M is such that 
M^I=e ?TL with 



q = ±1 if M has a real logarithm; 

(3.1.10) 

<j = ±2 if M has no real logarithm. 



The choice of sign of <j depends on the choice of the arrow of time. The linearized motion, even though 
it gives an idea of the nature of the true dynamics, cannot reproduce its fine structure. The time 
elimination makes the quadratic part of H integrable even though we have started from an a priori 
non-integrable one. The chaotic complexity is hidden in the higher-order terms and the aim of normal 
form theory is to retain only the main actors of this rich structure. 

3.2 Normal form for the quadratic part of the hamiltonian 

After having eliminated the time-dependence in the quadratic part of the hamiltonian, we are now 
ready to carry on the simplification of (p, q) by constructing a more suitable symplectic chart (p' , q'). 
The general classification of the normal forms of quadratic hamiltonians in any number of dimensions 
has been done a long time ago by 56] . In our case, since we are dealing with a bidimensional phase 



space, we will give a more hand-waving demonstration which however illustrates the main guidelines 
for deriving Meyer's classification in a pedestrian way. 

Let us characterize the corresponding symplectic transformation by a generating smooth func- 
tion S(p',q). The change of coordinates (p, q) i— * (p',q') is implicitly defined by 

fp=d q S(p',q); (3.2.1a) 

I q' = d p ,S{p\q) (3.2.1b) 
in the neighborhood of the origin, provided that 

d*, q S(0,0) (3.2.2) 

Since without loss of generality the transformation can be assumed to leave the origin invariant, 
S(p' , q) can be written as 

S(p',q) = \ 'Yp' 2 + 6p , q + ~r ] q 2 + O p , >q (3) (3.2.3) 

where (7,77) G 1R 2 and 6 £ IR\{0}, by virtue of condition ( |3.2.2| ). It is thus straightforward from 
equations ( |3.2.1j ) to obtain the linear part of the symplectic change of coordinates and deduce the new 
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expression of the quadratic hamiltonian 
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4o(c r? 2 + 2d 7/ + e) - 2 7 (d + c 77) + c 5 2 



/2 



def / 
= C 



--J-(cr/ 2 + 2dr/ + e) + d + cr7 



/ / 



dcf 



d' 



+ 



{erf + 2dr/ + e) 



<Z /2 +O pW (3) (3.2.4) 



def / 

= e 



Generically none of the coefficients (c, d, e) vanish. Nevertheless, we will also deal with situations 
where one real scalar like c, d, e or tr(M 2 ) — 2 can vanish a priori, since we are going to study not 
only individual hamiltonians but one-parameter unfoldings. Actually, we can choose (8, 7, 77) in order 
to cancel some of the coefficients (c',d',e'). For example, the cancellation of e' implies solving an 
algebraic equation of degree at most two in 77 whose discriminant is 



ce 



det(L) 



(3.2.5) 



The simplified form we want to obtain depends essentially on the nature of the eigenvalues of the 
monodromy matrix which governs the qualitative regime about the fixed point. 

Let us then briefly review the spectrum sp(M) of M. Since the monodromy matrix is symplectic, from 
condition ( |3.1,4 ) we have 



det(M) = 1 



(3.2.6) 



Moreover M is real and hence sp(M) is symmetric with respect to the real axis: v G sp(M) =>■ v* G 
sp(M). 

Let us now introduce the usual equivalence relation between matrices. Two matrices A and IB will 
be equivalent - we will denote A = B - if by definition they represent the same linear map but expressed 
in two different bases; in other words there exists an invertible matrix P such that A = P _1 BP. 
Without loss of generality P can be chosen with determinant one and thus symplectic. The Jordan 
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decomposition of M asserts that [51, chapter VI] 

where fj, £ R U LT(1)\{0} and ft € {0,1}. (3.2.7) 

It is possible to distinguish three different cases: 

(i) . sp(M) C IR\{— 1,0, 1} |tr(M)| > 2 and o is said to be an unstable fixed point; 

(ii) . sp(M) C {—1, 1} <=^- |tr(M)| = 2 and o is said to be a marginal fixed point 1 ; 
(in). sp(M) C U(l)\{-1, 1} -4=^ |tr(M)| < 2 and o is said to be a stable fixed point 2 . 

Generically, in the absence of any symmetry but if the dynamics is embedded in a one parameter 
family, the two eigenvalues of M may be degenerate (the marginal case). It would require at least 
two parameters to get M G {1, — 1} and, following our philosophy, we discard from the discussion all 
sets of codimension greater than two. In other words we will have \x = 1 in the generic marginal case, 



which can therefore be divided in two sub cases 



(ii) ' tr(M) = +2 and generically 
JJTj" tr(M) = — 2 and generically 



(1 ?); 
(-1-2) 




Figure 2: Geometrical representation of the locus of spM in the C plane. 



Moreover, case |^]] is usually also subdivided in the following way 
tr(M) > +2 and o is said to be a hyperbolic fixed point; 

[TJf ' tr(M) < —2 and o is said to be an inverse hyperbolic fixed point. 



(0 
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We now want to transpose this classification onto the spectrum sp(L) of L. Recall that the relation 
between L and M is given by ( |3 .1.1 0| ) . It can be easily verified that in cases [£)" and (11) " the mon- 
odromy matrix has no real logarithm and hence ? = ±2 while in all others cases - [[7J', [22)', \Hi\ - we 
can take <; = ±1. 



Since L is infinitesimal symplectic, property fl3,1.5|) implies that its spectrum sp(L) is symmetric 
with respect to zero: v G sp(L) —v G sp(L). Moreover, because L is real, sp(L) is symmetric with 
respect to the real axis: v G sp(L) v * G sp(L). Let 



L = 



A 

A -A, 



where 



AGRU(iR) and A G {0, 1} 



(3.2.8) 



be the Jordan decomposition of L. Then, relation (3.1.10) yields 



/ 



J2[2-|?I + (M -!)(/* + -)] 



\ 
1 



t — sinh(?r A) e 
* A 




(3.2.9) 



which implies that e 2rA = /i 2 with A G R U (ilR). Therefore, in the unstable case, we have A G K\{0} 
and then — det(L) > 0. In the stable case, we have Ar G i(R\7rZ) and then — det(L) < 0. In both 
cases ft = and A = 0. In the marginal case, we have M 2 = 1 + X where, generically, 



8),(-? a 8)} modulo = . 
Therefore, possibly changing the arrow of time, r i— > — t, we can choose 



(3.2.10) 




(3.2.11) 



then A = and det(L) = 0. Other choices, say L, of real infinitesimal symplectic matrices obeying 
relation ( 3.1.10| ) are possible and not necessarily with det(L) = 0. But following what we have done in 
subsection |3.l| , we can make a (linear and |?|r-periodic ) symplectic change of variables corresponding 



to F(t) = f e~'^e ^ in order to deal with an L whose determinant vanishes. Summarizing 
(i). unstable case A G IR\{0} and - det(L) > ; 
(it), marginal case A = and — det(L) = ; 

(Hi), stable case A G r _1 i (R\ttZ) and - det(L) < . 
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We can now go back to the reduction of the quadratic part of the hamiltonian ( fj.2.4 ). 



(i) Unstable case: We can choose at least one real value of r\ which makes e' vanish and, provided 
we take 7 = ±c 5 2 /[2^— det(L)], we can also cancel c'. At last, we get h^ 2 \p' , q') = ± \/ — det(L) p'q'. 
The global sign of the hamiltonian becomes irrelevant by a time reversal transformation and we finally 
find a symplectic chart in which the quadratic Hamiltonian can be written as 

hW(p,q) = \pq with A G R+\{0} (3.2.12) 

after having relabeled the new coordinates. h^ 2 \p,q) corresponds to an integrable dynamics whose 
trajectories lie on hyperbolae of center and whose asymptotes, the p and q axes, are the two stable 
and unstable manifolds respectively (see the right column of table ||). Recall (cf Eq.( 3.1.To| )) that in 
case (i) " the hamiltonian h^ 2 \p,q) describes the stroboscopic view of the real motion with a period 



twice the original physical period r. Technically such a prescription misses a — 1 factor in the equations 



of motion (3.1.1) and what effectively happens is that the system goes from one branch of a hyperbola 
to the other branch between t = and t = r. Then, between t = r and t = 2r, it goes back to the 
original branch. The 2r-stroboscopic view does not see this flipping. 



(ii) Marginal case: If we have d / and hence c / we can choose 77 = — d/c to get crj 2 + 
2dr/ + e = 0. Then both d' and e' vanish. Moreover, we take 5 = |tc| — 1//2 to obtain |c'| = 1/r. By 
changing possibly t to —t we get 

h {2) ( P ,q) = ^p 2 (3.2.13) 

after having relabeled the new coordinates. h^ 2 \p,q) corresponds to an integrable dynamics whose 
trajectories are straight lines parallel to the q axis (see the middle column of table ||). The distinction 
between cases JJTj' and JJT) n is analogous to that of the unstable cases. In (n) ', during one period r, 



the system flips back and forth between a couple of lines symmetric with respect to 0. 

(Hi) Stable case: We have cr/ 2 + 2dr/ + e 7^ for every real 7/ and then e' can never be cancelled. 



Nevertheless we can choose 7 = 5 2 (ct] + d)/(cr] 2 + 2d 77 + e) in order to have d' = 0. Besides, we 
can make e' = c' by choosing 5 = |(cr/ 2 + 2dr/ + e)) 1 / 2 det (L) -1 ^ 4 . If we define uj = |A| and possibly 
change the arrow of time we find the following normal form 

1 1 

h {2) (p,q) = -lo(p 2 + q 2 ) with w e - (R+\7t2) (3.2.14) 

2 T 

after having relabeled the new coordinates. h^(p,q) corresponds to an integrable dynamics whose 
trajectories are circles centered on and whose angular speed is to (see the right column of table |]). 

16 



in unstable case 



(ii) marginal case 



an 



(Hi) stable case 



M: monodromy 
matrix 



sp(M) C R\{-1,0,1} 



HI ?) 



tr(M) > 2 



tr(M) < -2 



JL/-1 o 



1 0\ 

1 -I) 



sp(M) c U(l)\{-1, 1} 



-2 < tr(M) < 2 



?€ {1,±2} 



? = 1 



kl = 2 



? = 1 



? = 1 



Ml f l = e ?TL 
sp(L)={A,-A} 



-det(L) > 



A e R+\{0} 



^ ~ 2t I 2 o ; 



det(L) < 



-iA G i (R+\ttZ) 



hV>(p,q) 



Xpq 



± v 2 

2tI ) 



flow around 
the origin 
( W = |A|) 












V 



P 




V 



Table 1: Classification of generic quadratic hamiltonians and normal forms in the neighborhood of an 
isolated fixed point o . 



Together with equations (|3.2.12|) , ( |3.2.13|) and it constitutes our first classification of normal 

forms. In the next subsections, we will refine it by including higher order terms in the hamiltonians 
and by considering one-parameter families of hamiltonians and normal forms. 



3.3 Birkhoff normal form for the full hamiltonian: unstable and stable cases 

If we write the equations of motion in complex coordinates 

z(p,q)=p + iq (3.3.1a) ( p(z , z) =^(z + z) (3.3.1c) 

z(p,q)=p-iq (3.3.1b) ( q(z, z) =-i(z - z) , (3.3.1d) 



then for the stable case (3.2.14) we get the differential equations z = Iujz and z = —iloz which correspond 
to the complexified hamiltonian h^ (z, z) = —iuj zz with (z, z) playing the role of complexified canonical 
variables. Moreover, from a function F(z',z) holomorphic with respect to its two complex variables 
in the neighborhood of the origin and satisfying <9 Z ' 5 5 -F(0, 0) ^ 0, one can construct a real symplectic 
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transformation of the form ( |3.2.1 ) for the generating function S(p',q) = — ip'q' + pq — $s [F(z', z)] V 



2 

where (p, q) and (z,z) are connected by the relations ( |3.3.1| ) and similarly for (p',q'), (z',z'). 
Then, we will study both unstable and stable cases from a hamiltonian of the form 

H(u,v;t) = Xuv + Y h " lc " 3 -""^e 1 " 1 ^* + O u<v (k + 1) (3.3.2) 

meZ 0,/3)eN 2 

a+/3 = k 

where we made use of the r-periodicity. 

From what precedes, if we are dealing with an unstable case then Ae1R + \{0}, (m, v) are the 
real coordinates (p, q) and the coefficients fi mQ(3 are complex numbers with a priori only one con- 
straint, ln* mal3 = h_ ma0 (H(p, q; t) must be real for real (p, q)). If the origin is stable, A G — i R + \{0} and 
(u, v) are the complexified coordinates (z, z) from which the real ones are obtained via equations fl3,3.1|) . 
In this case, H ma/3 are complex numbers such that W ma& = h_ m/3a (H(z, z; t) must be real for z = z*). 
In both cases, k is the least integer greater or equal to three such that h moi/3 / for some {a,f3) 
with a + f3 = k. 

We now want to simplify the expression of H by eliminating as many terms as possible in the 
series ( 3.3.2| ). The symplectic change of coordinates is given by the generating function 



S(u',v;t) = u'v + Y £ o ma0 u a v p e 1 ™^ + <v (k + l) (3.3.3) 
meZ («,/3)eN 2 

a+/3 = k 

where the transformation of degree less or equal to k — 1 corresponds to the identity since we do not 
want to modify the terms of the same order in H. a maj} are complex coefficients which we now fix. 
The symplectic transformation generated by the function ( |3.3.3| ) is 



' u(u',v';t) = u' + ^Y. E P^u^v'^e™ 2 ^ + O uW (k) (3.3.4a) 

meZ (a,/3)eN 2 
a+f3 = k 

v{u',v';t) = v'-^jY na^u^v'Ce™ 2 -? 1 + O uW (k) (3.3.4b) 

meZ (a,/3)6N 2 
a+f3 = k 
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dS 

and the hamiltonian H'(u' ,v';t) = H\u(u' ,v';t),v(u' ,v';t);t\ + -—\u',v(u',v';t);t] takes the form 

at 



H'(u',v';t) = Xu'v' + y~] ^ [h mQ(3 + (a(/3 - a) + im^)o- ma/3 



/a 1 8 im—t 



meZ (a,/3)eN 2 
a+8 = k 



+ O u /„,(fc+l) . (3.3.5) 

By taking cr mQ(3 so that the term between square brackets vanishes we can cancel all the terms in the 
double sum except the ones which correspond to m, a, 8 such that the so-called resonant condition is 
fulfilled 

2ir 

\{8-a) + im— = 0. (3.3.6) 
r 

Now the original hamiltonian may be simplified as follows. If we assume that all the non-resonant 
terms have been cancelled up to degree k — 1 in (p,q), we can get rid of those of degree k with the 
help of the generating function ( |3.3.3| ) provided we choose 

( v ma a = — ttt. r~~^ t f° r the- non-resonant m, a, B, 

Aid — a) + miZTr/T 

(3.3.7) 

I &m. a p = for the resonant m, a, 8. 

Since we can repeat this procedure order by order, we are left with a hamiltonian whose normal 
form is made up exclusively of resonant terms. This obstruction to the linearization of differential 
equations is well known from the works of Poincare, Dulac and later on Sternberg, Grobman and 
Hartman |2| among others. It is remarkable that the resonant relations which prevent the reduction 
via diffeomorphism are the same resonant relations as those we obtained above, even though they lie 
in a more restricted class of coordinate transformations (the symplectic ones). The convergence of the 
normalization procedure of a dynamical system is an old and delicate problem because of the possibly 
small denominators in expression ( |3.3.7 ). 



Let us now discuss the nature of the resonant terms depending on whether the origin is unstable 
or stable. 



JT) Unstable case: Since A G R + \{0} the resonant relation is fulfilled only if m = and a = 8. 
Thus, for n S IM\{0, 1} there always exists a real symplectic coordinate chart in which the hamiltonian 
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has the form 



H(p, q; t) =\pq + Y^ h k (pq) k + O p . q {2n + 1) 



(3.3.8) 



k=2 



for some 6 R. This normal form has been derived for the first time by || § III. 8] |25|] . n ^ 2 is 
a cut off integer whose value fixes the accuracy with which we want to reproduce the full dynamics 
and 2n will be called the order of the normal form. The limit n — > oo has been shown to be locally 
convergent p9], §4] and hence the dynamics is equivalent to an integrable one. 



1 



(m)| Stable case: We have A = ia> G — i([R + \7rZ). 
— r 



LOT 

2^ 



The hamiltonian can also be reduced to a Birkhoff normal form, which in real coordi- 



nates, is 



1 n 
H(p,q;t) = -uj{p 2 + q 2 ) + ^Q fe (/+9 2 ) fc + O p>q {2n + 1). 



(3.3.9) 



k=2 

for some £ (R and for every integer n ^ 2. In action- angle variables, 
1 



f / = | (p 2 + g 2 ) (3.3.10a) 



arctgl - 
\P 



(3.3.10b) 



' p = v 21 cos ( 
q = V2I sin ( 



(3.3.10c) 
(3.3.10d) 



and after having rescaled the coefficients, = f 2 k a.k, it takes the form, 



H(I,0;t) =loI + a k lk + 0/(2n + l) 

k=2 



(3.3.11) 



The problem of small denominators appears now clearly. The non-convergence of the normaliza- 
tion procedure as n — > 00 is due to the fact that the rational numbers Q are dense in R: cr mal3 in 
equation ( |3.3.7 ) can be arbitrarily large if we take (3 — a small enough. 

: Let t £ N\{0, 1, 2} be the smallest integer such that it exists r £ tN\{0} with 



7T 



r 2-7T 

u = tT 



(3.3.12) 



Then, the resonant terms correspond to (m, a, (3) £ Z x iNi 2 obeying (/3 — a)r = —£m. They are of 
degree s = a + (3 in (p, q). We can get rid of the time-dependence of the normal form if we describe 
the dynamics in a frame uniformly rotating around the origin at angular speed u. This corresponds 
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to the symplectic transformation 

( z i — ► z (z, z) = ze 



— it —\ — —iujt 

z [z,z) = ze 



generated by S(z',z) = z'ze lu)i , (3.3.13) 



leading to 

H'(z', z'; t) = Hfz'e-^, z'e iuJt ; t) + ^ (z',z'e iuJt ) (3.3.14) 

at v 



which has no longer any quadratic terms nor any time-dependence in all the resonant ones. 
The latter can be classified in three sets i?_ , Ro and R + defined by 

{(3 = a± k£ ( [3 = a 

R l for a = 0,1,2,... and k = 0,1,2,... . (3.3.15) 

s = 2a±n£^3 ^ s = 2a ^ 3 

Ro corresponds to the Birkhoff like terms. We obtain i?_ from R + exchanging (3 with a. 

These sets have a graphical interpretation in the (s,a)-plane (see figure |3| and Ref.|l3|]]). The 
resonant terms are associated with the points of integer coordinates s ^ a ^ and s ^ 3 which lie 
on straight lines of slope two, whose intercept are strictly positive, zero and strictly negative for 
Ro and R + respectively. 

For I ^ 3, in addition of the Birkhoff like terms, the resonant terms of lowest degree in (p, q) are 
(a = 0, (3 = t,m = —r) and (a = £, (3 = 0, m = r). Hence, dropping the primes of the rotating 
coordinates, we have 

H(z,z;t)= *k(zz) k + h* rto z e + h rm z e + O z>M (£+l), (3.3.16) 

ke N 

2^k^£/2 

for some E 1R and h rm £ C. In real coordinates (p, q), we get 

H( P ,q;t)= *k(p 2 + q 2 ) k + K M |K[(p + i^] + O m (£ + l). (3.3.17) 



feelM 

2<k<£/2 



In action-angle variables 



H(I,e;t)= a * J * + b^ /2 cos(^) + 0^1 + 1) (3.3.18) 

ke N 

2^k^£/2 

where = f 2 fc dfc, = f 2^/ 2 |h r( , |. In order to cancel the irrelevant phase 9q = f arg(h r£0 ), we have 
made a final rotation: 9 i— ► 6 — 9q. 
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a 



o: integer coordinates 



resonant terms 



Figure 3: Graphical representation of the resonant terms in a hamiltonian in the neighborhood of a 
stable fixed point whose linearized motion is a rotation of rational angular speed. 



3.4 Normal form for the full hamiltonian: marginal case 



We start from the real hamiltonian 



H(p,q;t) = —p 2 + E fWPVe im ^ + O p , q (k + l), 



(3.4.1) 



meZ (a,/3)eN 2 
ot+l3 = k 



k being a fixed integer greater or equal to three and h* ma/3 = H_ mQ(3 . We are looking for a suitable real 
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generating function of the form 



S(p',q;t)=p'q + £ £ a ma/3 p' a e™^* + Op/, ff (fc + 1) . 
meZ MeN 2 

a+/3 = k 



(3.4.2) 



Similarly to transformation (3.3.4) we get 



p(p',q';t) = p'+Y, E ^ 



cr mag p' a q" 3 - 1 e im r* + O pW {k) (3.4.3a) 



meZ (a,/3)GN 2 
Q +/3 = fc 



?(?/,(/;«) = d- E E «°™^P' a- V^e im ^ + O pU ,(fc) . (3.4.3b) 



meZ MeN 2 

o+/3 = fc 



For all integers m and (a, 0) such that a + (3 = k ^ 3, the coefficient of p' a q'@ e im ~* in the expansion 
of the hamiltonian H'(p, q';t) = H[p(p, q; t), q(p , q; t); t] + — [p', </; i); t] is 



27T 

h ma/ 3 + im & ma p — 2(a + 1) a m r a+1 - ) m_ 1 - ) 

T 

■ 271 

r 



Va^0,V/3>l such that q + (3 = k ; 



(3.4.4) 



for a = and (3 = k. 
If we want to cancel them for fixed and m, we must solve the matrix equation 



im 



2- 



im^ 




-4 







im 



2?r 



-2(a + l) 



im 



2tt 





-2fc 



im 



2tT 



(J„ 



(J, 



ml(fc-l) 



(7, 



<7„ 



h, 



L ml(fc-1) 



(3.4.6) 



which we denote by 



Aq a = H 



(3.4.7) 
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Since det^o = (im^) fc+1 , the above equation has always a solution for every m 7^ and h, thus 
all the time-dependent terms can be cancelled. When m = 0, the kernel of Aq is one-dimensional 
and generated by Ui = f (1,0, ... , 0) and the image Im Aq is a hyper plane of equation Ufc+i-h, = 
with Ufc + i = f (0, . . . , 0, 1) and where we have introduced the canonical scalar product. Let us now 



decompose the right hand side of Q3.4.6 ) in the following manner 

h = h|| + h ± (3.4.8) 

where fly £ ImAo and h. ± .Im^4o = 0. We can choose a = A^ 1 ^, and the only coefficients which 
cannot be eliminated correspond to the elements orthogonal to Im Aq , that is those which are parallel 
to Ufe+l 

h-oofc = h 01 ( fc _ 1 j = • • • = \ioa(k-a) = ' ' ' = h-o^.^o = and H 0fc0 ^ . (3.4.9) 

The latter are therefore the resonant terms and the form of the hamiltonian is, after relabeling the 
coordinates 



1 n 

H(p,q;t) = —p 2 + ]Th fc g fc + M (n + l), (3.4.10) 
for some h.^ £ [R and for any integer n ^ 3. 



2t 

k=3 



4 Normal forms for one-dimensional parameter families of hamilto- 
nians 

Following what we have planned in section ||, we will now consider that the hamiltonian H (p, q; t) is 
embedded in a smooth family depending on one real parameter, say s. A typical example of the control 
parameter is the energy of a non-reduced two degrees of freedom system. Without loss of generality, 



we assume that the origin o is an isolated equilibrium point of H(p,q;t;e), obeying equations ( 2.0.4 ) 
at e = 0. We can follow step by step the discussion of the previous section for H(p, q;t;e = 0) and 
classify its generic normal form around 0. Now the question is to study what happens to the dynamics 
if we vary e in a neighborhood of zero. We must first wonder about the behavior of the fixed point. 
By virtue of the implicit function theorem applied to the equations 
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Vt E R, 



<9# 

dH 
dq 



(p,q;t;e) = 
(p,q;t;e) = 



(4.0.11a) 



(4.0.11b) 



in a neighborhood of (p, q; e) = (0, 0; 0), we can follow as a function of e the isolated solutions if 

/ dl v H{w; t; e) d*H(w; t; e) 



+ det (Hess [#(«;; i;e)]). 



=(°),*=o 



det 



)2 U („.,.+. r-\ »2 



\dq )P H(w; t; e) d q z )q H(w; t; e) 



. (4.0.12) 

An equivalent condition may be written for the monodromy matrix M. If w i— > $(u;;e) is the (local) 
Poincare map, the implicit function theorem applied to 



<£(u>; e) = w 

states that we can follow smoothly the fixed point if 

^ det[d w ^(w;e)-t), =(0)te=0 = det [M(e) - l] |e=0 



tr [M(e 



I | E =o • 



(4.0.13) 



(4.0.14) 



The last equality comes from property ( p.2.6|) . From table |l] we see that the latter condition is fulfilled 
in both stable and unstable cases as well as in case Kii)f , because there exists a neighborhood J? 
with 0G/cK and a smooth family of points p(e) such that p(0) = 0, and p(e) is a fixed point 
of $(-;e) for all e in J^. If we now make an appropriate smooth e-dependent change of coordi- 
nates in order to translate back p(e) to o, then we deal with a hamiltonian family H(w;t;e) which 
admits o as equilibrium point for all £ in /. Elimination of the time dependence of the quadratic 



part of the hamiltonian around o can be done as described in subsection 3.1 in the whole . The 
r or 2r-periodic moving frame in which we follow the dynamics is smooth with respect to e. At this 



stage, we must distinguish cases (£) and {in) from case (ii)" . On the one hand, in the former two cases 
the stability of the fixed point will not change if ^ is small enough because from a topological point 
of view property tr[M(e)l ^2 defines an open set for e. On the other hand, case Ju) n corresponds 
to tr[M(e = 0)] =—2 which breaks down as soon as e^0. More precisely, generically tr[M(e)] will 
cross transversally the value —2 and o will change from stable to unstable being marginal only at e = 0. 
We will study first the stable and unstable unfoldings, then the marginal unfolding will be considered 



and we will finish with the study of case (ii) ' where even the presence of a fixed point cannot be 
guaranteed as soon as e ^ 0. When embedded in a one-parameter unfolding, case and |(m)|" are 



respectively called extremal and transitional by [34, resp. definition 1.9 and definition 1.12]. 
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4.1 Normal forms for unstable and stable unfoldings 

From what precedes, we can deduce that if o is an unstable (resp. stable) fixed point of H(w;t;e = 0) 
it will remain a fixed point of H(w; t; e) for all e in a small enough neighborhood J' of zero. Indeed, 



in both cases all the symplectic changes of coordinates we have made in subsection |3,2| and |3.3| vary 
smoothly with e. 

In the unstable case, we are led to the Birkhoff normal form (3.3.8) with now A and {^fcj^g^ n} 
being smooth real functions of e 

n 

H u (p, q; t; e) = X(e) pq + Y, M^) (rf + O m (2n + 1) . (4.1.1) 

k=2 

IT 

In the stable case, we cannot prevent uj from crossing — Q\tNl infinitely many times. Let us then 

r 

assume that 

a;( £ = 0)=u,o^ f ~- (4.1.2) 

£ T 

with r G N\{0} and £ G N\{0, 1,2} coprimes. Choose the neighborhood / D so that it contains 
no rational ^uj{e) whose denominator is lower than I. Moreover since generically ^| £ _ 7^0, we can 
always choose the control parameter in order to satisfy uj(e) — loq = e. Thus, in a rotating frame 
defined by ( p.3.13j ) with uj = ujq we get, from expressions ( 3.3. 17[ ) and ( 3.3.18| ), 



H e (p,q;t;e) = \e{p 2 +q 2 )+ £ a fc (e) (p 2 + q 2 ) k + \h m (e)\ R[(p + iq) £ ] + O m (£ + 1) 

(4.1.3a) 



2^fc^/2 



and 



H e (I,8;t;e) = el + a k (e) I k + b £ (e) I l/2 cos(^) + 0^(1 + 1) (4.1.3b) 

ke N 

in momentum-position and action-angle variables, respectively. 



4.2 Normal form for the transitional fixed point 

Since the stability of the fixed point will change when crossing e = 0, derivation of the normal form 
of the unfolding is not as straightforward as in the previous subsection. Let us then start from the 
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following hamiltonian 



H(p, q; t;e) = - c(e) p 2 + d(e)pq + - e(e) g 2 + O m (3) . (4.2.1) 



= f h {2 \p,q;e) 
c, d and e are smooth functions of e such that 



c(0) = - , d(0) = , e(0) = . (4.2.2) 
r 



In Eq.(4.2.1) there are no linear terms in (p,q) because we have already shown that the origin o is 



fixed in a suitable symplectic chart for all e. Recall that the hamiltonian ( [4.2.1 ) corresponds to a 
2r-periodic stroboscopic view of the true dynamics. 



We will proceed as in subsection 3.2. We look for a family of new symplectic charts obtained via 



a family of generating functions of the form 

S(P, T,e) = \ 7 (e) P 2 + 6(e) p 'g + \^) q 2 + <V, , (3) (4.2.3) 

where (5,j,r]) are now smooth real functions of e such that 5 does not vanish in ^ . We are led to 
an expression identical to ( ggg ) for h^(p',g'). The length of J can possibly be reduced in order 
to have, Ve G J, c(e) > 0. Moreover since generically rf ^ d d J ce ^ ^ £ _ Q 7^0, we can reparametrize the 
unfolding by e' = — r(d 2 — ce). We can then choose 



We el, 7(e) = 0; 6(e) = l/y/rc(e); 17(e) = -d/c , (4.2.4) 

which implies that c' = 1/t, d' = 0, and e' = e. 

Dropping all the primes, the normal form of the quadratic terms of the unfolding of a transitional 
point is 

hW(p,g;e) = ^p 2 + eq 2 . (4.2.5) 
It 



We can now reproduce step by step the reasoning of subsection p.4| . Starting from 

H(p,q;t;e) = ^-p 2 + eq 2 + E p a g P e™^ + O p , q (k + 1), (4.2.6) 

meZ («,/3)eN 2 

we construct a generating function ( |3.4.2| ) with the cj's being smooth functions of e. We are led to 
a matrix equation A(e) a(e) = h(e) with ^4(e = 0) = Aq [see ( 3.4.6j )1- Possibly reducing the length 
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of J^, we can assume that Ve E ^,det^4(e) 7^ when m/0, and then the time dependence can be 
cancelled uniformly in e. In fact, we can also eliminate all the monomials of order greater than three 
except those terms corresponding to ( [3.4, 9 ) since we cannot cancel them uniformly in e. We thus 



obtain the full normal form for the unfolding in the transitional case 

n 

H(p,q;t;e) = —p 2 + eq 2 + Y,^k(e)q k + O m (n + l), (4.2.7) 



2r- 

k=3 



where {hfc}fce{3,...,n} are some real functions of e. 



There are certain constraints on the H's due to the fact that the hamiltonian (4.2.7) describes the 



dynamics in a (— 2"r)-periodic stroboscopic view. After one period r, the monodromy matrix verifies 

l(e = 0)= ( "| I . (4.2.8) 



It can be shown using the Lie transformation [37, §VII. E. 5] that we must have H(p, —q;t;e 



0) = H(p,q;t;e = 0) + VA (n + 1) which implies \/ k e N\{0}, h.2fc+i(0) = 0. We are going to prove 
the same result for the first terms of the normal form (up to n = 6) but without the help of the Lie 
transformation theory. We will show in this subsection that 113(0) = 0. We will later on show in 
subsection that ^ri e=0 = and h,5(0) = (see equations ( |6.3.10| )). 

Let us fix e = 0, define V3(q) = h^(0)q 3 and anticipate a little bit on the calculus of Poincare 
stroboscopic maps of section |(| We apply formula ( 6.0.2j ) in order to calculate the Poincare map at 
time —It. The Poisson brackets are given by 

\q,H\ = -p + O p ,g(3) ; (4.2.9) 
r 

{p,H} = -v' 3 (q) + O p ,,(3) ; (4.2.10) 

Up,H},H} = -^pvUq) + Op,,(3); (4.2.11) 

{{{p,H\,H},H} = -lp 2 4"(g) + O p , q (3) ; (4.2.12) 

{...{p,H\,...,H} G Pi9 (3) . (4.2.13) 

V v ' 

n times, n^4 

Then it is straightforward to get W2 = w(—2t) up to order two in wq = w(0) 
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p 2 =p + 2Tv' 3 (q ) - 2Tp v 3 (q ) + - TPo^'fao) + Po ,go( 3 ) I 

4 2 

k 92 = go - 2p - 2t*4(<7o) - g Tpo^'Cgo) - , Tpg^"(g ) + Po ,go( 3 ) ; 



(4.2.14a) 



(4.2.14b) 



3 JV ' 3 

Now we use the fact that this map is the square of another map whose linear part is known (see 
equation Q4.2.8Q ) 



' P\ = -Po + /2(Po,9o) +O p(hqo (3) ; 



(4.2.15a) 



Ml = -Qo +P0 + 92(po,Qo) + O po>qo (3i) ; (4.2.15b) 

(/2 and 52 are two quadratic forms in w). After having iterated the last map and identified the result 
with (14.2.14| ) we get 

2ru£(?) - 2rpt4'(g) + ^rp 2 <(<?) = f2(-p,P~q) - f 2 (p,q) (4.2.16) 

for all (p, q) in a neighborhood of the origin. In particular on the g-axis relation ( [4.2. 16[ ) implies h,3(0) = 
0. We have then shown that the unfolding of a transitional point has the form 

1 " 
H(p, q;t;e) = —p 2 + sq 2 + eh 3 {e)q 3 + h*(e) ^ + M (n + 1), (4.2.17) 

were h^e) = f h.3(e)/e is a smooth real function. 



4.3 Normal form for an extremal fixed point 



In this case condition ( 4.0.14 ) is not fulfilled. Generically we cannot guaranty that we can follow 
smoothly the fixed point as soon as e / and then we cannot even make use of the arguments of 
subsection |3. 1| to remove smoothly the time-dependence of the quadratic part of H (p, q'; t; e) uniformly 
with respect to e in a neighborhood of zero. We must therefore reconsider our approach from its very 
beginning. 



4.3.1 Uniform elimination of the time-dependence in the quadratic part of the unfolding 

The most general form for a one-parameter linearized r-periodic differential equation is 



w(t; e) = A(t; e) w(t; e) + a(t; e) 



(4.3.1) 
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where A(t; e) (resp. a(i; e) ) is a real one-parameter family of T-periodic 2x2 matrices (resp. 2- vectors). 
If U(t;e) is the operator defined by equations ( 3.1.6| ) and ( 3,1.7| ), then the extremal case corresponds 



precisely to U(t;e = 0) = M(e = 0) = (\ 5). From now on and until the end of this subsection 
we should keep in mind that all the quantities we are dealing with depend smoothly on e but for 
simplicity we will not write down explicitly this dependence. Let L (resp. b) be any smooth family of 
time-independent infinitesimal symplectic matrices (resp. 2-vectors). In the new symplectic coordinate 



system defined by w(t) = e*^ [U(i)l 1 w(t) + b(t), equation ( 4.3.1 ) becomes 



w(t) =hw(t) +b-hb + e [U(t)] 1 a(t) . (4.3.2) 

* v ' 

def „ 

= a 

There will be no explicit time-dependence if 

b-~Lb + a(t) = c (4.3.3) 

where now c is a (family of) constant 2-vector(s). As we have seen in subsection [O], the change of 
variables has to be real and r-periodic. Actually it is an extremal case if we choose a real r-periodic b 
and L = = ( J q ) • Let us then expand a(t) and b(t) in Fourier series 

a(t) = Sne in ^* and b(t) = ^ 6„e in ^* where Vn G Z, (a n ,b n ) £ C 2 . (4.3.4) 

neZ neZ 



Then equation (|4,3.3j) leads to 



VnGZ, (in^l-L) b n = 5 ( n K ^c- a n (4.3.5) 



T 



(<5 (Jf) is the Kronecker symbol). Since the eigenvalues of L are both small for e near zero, the matrix 



equation (4.3.5) can be solved for every n E Z\{0}. It only remains L&o = ®o ~ c which can also be 
solved provided that So — c belongs to ImL = 1R(q) +O e (l). If ao = {al) we can therefore eliminate the 
time dependence of ( 4.3. 1| ) uniformly in e by choosing c = (£ ) + O s (l) and b = L-l(a - c) + O e (l). 



Moreover, since L is real we have Vug Z, 6* = 6_ n which ensures that b(t) is real. 

We have then shown that even in the one-parameter unfolding of an extremal point, we can choose 
a suitable local symplectic coordinate system smooth in e in which the dynamics corresponds to a 
hamiltonian whose terms of degree one and two in (p, q) are constant. 
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4.3.2 Normal form for the quadratic unfolding 

From what precedes, let us start from the most general form for a time-independent hamiltonian 

(4.3.6) 



H(p, q; t; e) = a(e) p + b(e) q + - c(e) p 2 + d(e) pq + - e(e) q 2 + O m (3) 



M^( 2 )(p,g;e) 

with (a, b, c, d, e) being smooth functions of e such that 



a(0) = ; b(0) = ; c(0) 



1 



d(0) = ; e(0) = 



(4.3.7) 



If .y is small enough we have Ve G J^, c > 0. Then, if we make the smooth translation p = p-\ — 

c 



and if we relabel the coefficients of , we can write 



hF*(p,q)e) = b(s) g + - c(e) p 2 + d(e) pq + - e(e) q 2 . 



(4.3.8) 



As was done in subsection |4.2| , the symplectic change of variables corresponding to the generating 
function ( [4.2.3 ) gives 



h {2 \p',q';e) = \hq' + \ 



^{crf + 2d ?? + e) - 2 7 (d + erf) + c<5 2 



12 



(4.3.9) 



+ 



(cr? 2 + 2dr/ + e) + d + cr] 



pq + 



{erf + 2dr/ + e) 



9 ,2 + P ' l9 '(3) • 



Let us make the same choice as in (4.2.4). Aside from Eq. fld.3.7 ), generically we have no addi- 



tional relations between (a, b, c, d, e), in particular 



de [e=0 



^ and without loss of generality we can 



reparametrize the unfolding to e' = b-^/rc. Dropping all the primes, the normal form of the quadratic 
terms of the unfolding of an extremal point is 

h {2) (p, q;e) = ^p 2 + eq + eh 2 (e) q 2 , (4.3.10) 

where h.2 is a smooth generic real function of e. Getting the normal form for higher orders can be 
done uniformly with respect to e exactly in the same way as described in subsection 4.2. A full normal 
form for the unfolding in the extremal case can then be written 

1 n 
H{p,q;t-e) = —p 2 + eq + e\i 2 {e)q 2 + Y,^k{e)q k + O p>q {n+l), (4.3.11) 



fc=3 



where {bfc}fce{3,..., n } are some generic real functions of e. Through the smooth translation q' = q + |pj^j 
we can moreover set ^(e = 0) = 0. 
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5 Geometrical structure of the trajectories 



We have now all the necessary elements to describe the unfolding of the dynamics for each of the 
above mentioned hamiltonian normal forms. In particular, staying in a neighborhood of the origin o 
we search for a characterization, as a function of e, of the typical phase-space scales of the classical 
structures involved in each bifurcation. We also want to explicitly locate the real as well as the complex 
fixed points involved. The following will be valid as long as e remains small compared to all algebraic 
expressions that can be made from the coefficients appearing in the normal form and having the same 
dimensions. 

Recall that the dynamics generated by the normal forms coincides with the true one if we look at 
the trajectories stroboscopically with a period that depends on the nature of the equilibrium point 
at e = (see table [l]). We must also keep in mind that, in the stable case, the hamiltonian ( 4.1.3| ) 
describe the dynamics in a rotating frame (see equations ( 3.3. 13| )). 



5.1 Unstable Point 



The unfolding of an unstable point is given by hamiltonian ( 4.1.1 ). For each e in J?, the origin o is 



the only fixed point and it remains unstable for J 1 small enough. 



5.2 Extremal bifurcation 



The normal form (4.3.11) can be rewritten as follows (recall that we can take h.2(£ = 0) = 0) 



H x {p,q\t;e) = ^p 2 + eq + ^aq 3 + l -bq A + O p , 9 ,^(5) (5.2.1) 

where a, b are generic real numbers and e is such that |e| <C \a 3 /b 2 \. In figure || we have plotted the 
potential f ex tr(<?) = £q + § aq 3 whose behavior governs the dynamics when q ~ 0. The fixed points 
of Hi(p, q; t; e) are given by p = and q 2 = —e/a. 



e/a < There are two fixed points which are real 

/ 



def 
P± = 



p = 

def 



\ q = q± = ±vWo| 



+ 0^(2) . (5.2.2) 



The equations of motion are 
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V 



-e - aq 2 + 0^,^(3) ; 



,g= p + O p ,5,^(3) 

Their linearization around p± gives 





(5.2.3a) 
(5.2.3b) 



- 2aq± \ ( u 
1/r )\v 



(5.2.4) 



Applying the results of subsection 3.2 we deduce from the discussion of subsection 3.2 that psgna is 



stable with a local angular speed given by uj = A/2/r|ea| 1//4 . p- sg na is an unstable point whose 
stable and unstable manifolds have respectively the following tangents: u = —y / 2/r\ea\ 1 ^v 
and u = y / 2/r|ea| 1 / 4 f . The equation of the separatrix is 



^p 2 + eq + ^aq s = Hi(0,q- sgna ;t;e) 



3 

£ 2 



(5.2.5) 



It is the boundary of an area whose height and width vary with \e\ as \e\ 1 ^ 2 and |£| 3//4 , respec- 
tively (see figure |5|) . 



e/a = By construction, we have one and only one fixed point: the origin o. The two fixed 



points p± which are present when e/a ^ now coalesce. The linearized flow has been studied in 
section 3.2 and looks like the marginal case in table |]. If we take into account the first non-linear 
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q 




Figure 5: Flow near the bifurcation of an extremal fixed point when e/a < 0. 

q 




v 



Figure 6: Flow near the bifurcation of an extremal fixed point at e = 0. 



terms we get two families of smooth trajectories separated by a curve given by 



J-p2 + I a g3 = 0, (5.2.6) 



which has a vertical cusp at o where the stable and unstable manifolds meet (see figure 

34 



e/a > There are two fixed points which are now complex 



P± 



dcf 



p = 



(5.2.7) 



q = i q± = ±iy/\e/a\ 

In the Poincare surface of section, the flow has no fixed point and the linearized flow is given by 
equations ( |5.2.4j ) (see figure ^). 




V 



Figure 7: Flow near the bifurcation of an extremal fixed point when e/a > 0. 



5.3 Transitional point 



The normal form ( 4.2.17 ) can be reordered as follows 



1 



H 2 (P, q; t; e) = — p 2 + eg 2 + - a + - e b q A + - eg 5 + h4(0)e g A + h 6 (0)q b + ^_(7) 



(5.3.1) 



where a = 4h 4 (0), b = 3h 3 (0) = 3h' 3 (0), and c = 5h 5 (0). We are in a regime where |e| <C \a/b \ 
and |e| <C |c/fo 3 |. In figure |8| we have plotted the potential Vt Tan (q) = £Q 2 + j^q 4 whose behavior 
governs the dynamics. 

The fixed points of H(p, q; t; e) are the origin o and the two points defined by p ~ 0, q 2 ~ —2e/a. 



dcf 



dcf 



e/a < Aside from the origin, there are two other real fixed points 



P± 



dcf 



p = 



Vg = & d = f ±v / ^R7 



(5.3.2) 
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Figure 8: Potential wtran(?)- 



The equations of motion are 



p = -2s q - a q 3 + q> ^ (4) ; 



1 «=r P + °V^uM A) - 



(5.3.3a) 
(5.3.3b) 



Their linearization around o gives 





1/r 



- 2s 




(5.3.4) 



Therefore o is stable if a < and the angular speed around it in this case is to = 2y / 2|e|/r ■. 



If a > 0, o is unstable and the stable and unstable manifolds have tangents u = ^fy/2\e\/r v. 



The linearization of equations (|5.3.3[) about p± gives 



As\fu 
1/rO )\v 



(5.3.5) 



Then p± have a stability opposite to o: if a < 0, p± are both unstable and in a chart centered 
on them, the stable and unstable manifolds have tangents u = ±2^\s\t v. If a > 0, p± are both 
stable and the angular speed around them is uj = 2y / \s\r. 
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The equation of the separatrix is 



1 



1 



2T P 2 +eq z + -aq* + O VWl (5)=H 2 (0,0;t;£)=0 



if a > 



(5.3.6) 

, A aq i + O VWi (5) = H 2 (0,q ± ;t ] e) = -e 2 if a < . 

In both cases, the separatrix is the boundary of an area whose height and width varies with |e| 
as |e| 1//2 and |e|, respectively (see figure |9|). 

Q q 




P 




P 



(a) a < (b) a > 

Figure 9: Flow near the bifurcation of a transitional fixed point when e/a < 0. 



e/a = By construction, the origin is the only fixed point. The two fixed points p± which are 



present when e/a ^ have merged at the origin. The linearized flow has been studied in section 
3.2 and looks like the marginal case in table [j]. If we take into account the non-linear terms we 
get families of smooth trajectories separated by the curve ^p 2 + \aq A + 0^1(5) represented by 
a couple of tangent parabolae if a < and a single point if a > (see figure |To| ). 



e/a > o is the only real fixed point and there are two fixed points in the complex plane 

\ 



P± 



dcf 



p = 



\ Q = k± = ±iy/2\e/a\ J 



+ <W 2 ) • 



(5.3.7) 



The equations of motion and their linearization at o are given by Eqs. (5.3.3) and ( [5.3. 4 ). 
respectively. The origin changes its stability when e crosses zero, o is stable if a > and the 
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Q 




f V 



Figure 10: Flow near the bifurcation of a transitional fixed point at e = 0. 



angular speed around it is u = y / 2\e\/r. If a < 0, o is unstable and the stable and unstable 
manifolds have tangents p = =Fy / 2|e|/r q (see figure 11). 

q 




p 




(a) a < (b) a > 

Figure 11: Flow near the bifurcation of a transitional fixed point when e/a < 0. 



As we have said (see table |l]), the — 2r-periodic stroboscopic view of these flows and the doubly 
iterated Poincare map coincide. We will see, from the study of the generating functions of next section, 
that in order to describe the true dynamics after one-step iteration of the Poincare map we must have 
a global symmetry with respect to o. That means that p± are non- longer fixed points when e ^ but 
rather belong to a periodic orbit of period 2r. In one period, they exchange each other and all the 
points in their neighborhood jump across o. 
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5.4 Stable point with u = 2m /(It) and £ ^ 3 
5.4.1 £ = 3 

For I = 3, the hamiltonian ( |4.1.3| ) takes the form 



H 3 (p,q;t;e) = \e{p 2 + q 2 ) + - b(p 3 - 3pq 2 ) + O p , ?)£ (4) 



and 



1 



H 3 (I,6;t;e) = el + ^ b (2/) 3 / 2 cos(3 9) + 0^ £ (4) 

for some generic real number b. The fixed points of ( [5.4. 1| ) are real for all e. Aside from the origin we 
have 



(5.4.1a) 



(5.4.1b) 



Po 



dcf 



p = —e/b 



( 



dcf 



p = e/{2b) 



\ q ± d ^±V3e/(2b) 



+ £ {2) . (5.4.2) 



+ O e (2) and p ± 

\q = 

These points lie near the circle centered at o with radius \e/b\. From expression ( 5.4.1b| ) we see that 
the flow is - to first order - invariant by rotation of around o. The equations of motion are 



p=-eq + 2bpq + Pi g j£ (3) ; 



(5.4.3a) 



q = ep + b(p 2 -q 2 ) + O m , £ (3). (5.4.3b) 
By construction (see subsection |4.1| ) about o we have a rotation of angular speed e. The linearization 



of ( 5.4.3 ) around po gives 





-£ 



3e\ (u 
J U 



which shows that po an d p± are unstable. The separatrix has the equation 



= 6 



2b 



1 / E 



[ P+ b 



1 / e 



(5.4.4) 



(5.4.5) 



which, to first order, is made by the three straight lines defining the equilateral triangle (po, p+, P-) (see 
figure |l2|-a). At e = 0, all three points collapse at o (see figure |l2|-b). 

The flows are seen in a frame rotating at angular speed wo = 27rr/(3r) (cf (|3 .3.1 3| ) ) . Thus, in the 
fixed frame, {po,p+,p-} are no longer fixed points but rotate with angular speed wo- They belong to 
a periodic orbit of period 3r. 
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q 




sgn(f)p 




p 



(a) e # (b) e = 

Figure 12: Fiow near the bifurcation of a stable fixed point with t = 3 



5.4.2 £ = 4 



For £ = 4, the hamiltonian ( |4.1.3| ) takes the form 



H A (p,q;t;e) = - e (p 2 + q 2 ) + - a(p 4 + q 4 ) + - bp 2 q 2 + O p>WiiI (5) 



(5.4.6a) 



and 



H A (I,e-t;e) = el + -(3a + b)I z + -(a - b)I l cos(4 6) + 0(^^(5) 



(5.4.6b) 



for some generic real numbers a and 6. Because we can change the global sign of H by inverting 
the arrow of time and the sign of e, we can consider a > without loss of generality. Moreover, by 
possibly make a global rotation of 7r/4, we can assume a — b > 0. From expression ( 5.4.6b| ), we see 
that the flow is - to first order - invariant by rotation of around o. The fixed points are solutions 
of the algebraic system 



or, in action-angle variables 



ep + ap 3 + bpq 2 + P:<1 ^{A) = 
eq + ap 3 + bp 2 q+ O^^A) = 



(5.4.7a) 
(5.4.7b) 
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'sin(40) + O v7 ^(3) = O; 

, e + - [3a + b + (a - b) cos(40)] J + ^(3) = . 



(5.4.8a) 
(5.4.8b) 



Since generically a / 6, we must distinguish the two cases a + b ^ 0. 



a) a + 5 > 



e/a < There are nine fixed points which are all real. Aside from the center we have 



± def 



p = ±y / -e/a 

\ q = o 



P± = 



P = ±^-e/(a + b) 



\ 



modulo Pj?)V ^|(2). 



o = ±V-e/(a + 6) / 



i fa def 

and p± = 



( p = \ 



\q = ±y / -e/a J 



(5.4.9) 




P 



Figure 13: Flow near the bifurcation of a stable fixed point of £ = 4 when a + b > and e/a < 0. 



Linearizing the equations of motion we find that the four points {p^} are stable while {pijp 1 * 1 , } 



41 



are unstable. The separatrix between them is given by 



i(p 2 + Q 2 ) ~ ^2a{a-b)pq + ej \a{p 2 + q 2 ) - ^2a{a-b)pq + ej + 0^,^(3) = . 

(5.4.10) 

To first order the separatrix is made of two ellipses whose larger axes are oriented along 9 = 
±7r/4. Their intersection defines four satellite islands of typical width proportional to yje| (see 
figure liH) . 



e = is the only fixed point and it is stable with a small positive angular rotation around it. 



The flow looks like the stable one in table [I]. 



e/a > o is the only real fixed point and it is stable with angular rotation e around it. There 

/ p = 



are also eight fixed points with complex coordinates 



/ 



± def 



p = ±iy e/a 



± def 



/ p = ±iy/e/(a + b) 



\q = 



and p± 



def 



\q = ±Ve/(a + b) 



i q = ±.\yje/c 



(5.4.11) 



modulo O p q ^^{2). The flow looks like the stable one in table |. 



(3) a + b < 



e/a < 



There are five real fixed points 



o, p 



± def 



p = ±sj-e/t 
\q = 



( 



, and p± 



def 



p = 



\ 



\ q = ±^-e/a J 



(5.4.12) 



modulo O p q y ^(2), and four fixed points with complex coordinates 



f p = ±iy / -e/(a + b) 
\q = ±iy / -e/(a + b) 



modulo P)9iv ^|(2). 



(5.4.13) 
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Q 




V 



Figure 14: Flow near the bifurcation of a stable fixed point when a + b < and e/a < 0. 



Linearizing the equations of motion we see that the four points {pijp 1 * 1 , } are all unstable. The 
separatrix joining them is given by 



a(p 2 + q 2 ) ~ V 2a ( a ~b)pq-e a{p 2 + q 2 ) + ^2a{a -b)pq-e + Op, giV ^(3) = . 

(5.4.14) 

To first order, it is locally made of two hyperbolae of principal axes 6 = ±7r/4 (see figure |j~4| ). 



e = The nine fixed points have collapsed into a single point o which is unstable (see figure |lq) . 



e/a > 



In addition to the stable fixed point at the origin, there are four real fixed points 



± def 
P± = 



p = ± v / e/(a + b) 



(5.4.15) 



q = ±^/e/{a + b) 

modulo Pi(?jv ^(2). They are all unstable. To first order, the separatrix between them is 
given by two hyperbolae like in Eq. (5.4.14) but now the principal axes are oriented along 9 = 
and 9 = ir/2. Besides them, there are four fixed complex points 
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Q 

L 




Figure 15: Flow near the bifurcation of a stable fixed point of I = 4 when a + b < and e = 0. 



o, p 



± dcf 



V 9 = 



and p± 



(U'f 



p = 



\q = ii^-e/a y 



See figure |T(]. 




(5.4.16) 



Figure 16: Flow near the bifurcation of a stable fixed point when —b > a and e/a > 0. 

In all cases described in a) and (3) the reference frame rotates at angular speed ojq = 7rr/(2r) (cf 
Eq. 3.3.13 ). In the fixed frame, the points p rotate uniformly around o with speed Uq. {p^} belong to 



the same periodic orbit of period 4r. 



5.4.3 £>5 



We will work in action-angle variables with the normal form ( 4.1.3b ) . The origin remains stable across 
the bifurcation but the sense of rotation of the flow around it is inverted. The other fixed points are 
solutions of 
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- £\)il e / 2 sm{£9) + 0^^(1+1) = 0; 

e + Y, ka k (e)I k ~ l + \i b e (e) f' 2 ' 1 cos(^) + 0^^{i-l) 
ke N 

To first order they correspond to: 9 k = f k— for k £ {0, ...,2^— 1} and / 
where a 2 = a2(e = 0). 



(5.4.17a) 
. (5.4.17b) 



e/cb2 < In addition to the origin, there are 2£ real fixed points 
A J 1 = -e/(2a 2 ) \ 



Pk 



\ k d = kir/i ) 



+ 0^|(3) with k £ {0, ...,2^ — 1} . (5.4.18) 
To first order, they lie on the circle with center o and radius \J—e/a2- The linearization of the 




Figure 17: Flow near the bifurcation of a stable fixed point for £ = 5 when eja^ < 0. 
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equations of motion about p^ is 

' r 



{-l) k e 2 b\e/2a 2 \ i/2 \fl 
2a 2 



(5.4.19) 



which shows that the points {pk}ke{o,...,2£-l} are alternately stable and unstable. Without loss 
of generality, we can assume that po is stable. The separatrix has the following equation 



el+ a k (e)I k + b £ (e)I e / 2 cos(£0) = H i (p 1 ;t;e)+ + 



(5.4.20) 



where 



dcf 



if £ = 3 ; 
if £ + 3 . 



(5.4.21) 



It is a bounded set ^-periodic in 6. To evaluate the typical transversal size Sr of the island 
around the stable satellites we should solve to first order equation ( 5.4.20| ) for 9 = with 
/ = /(pi) + SI. We obtain 51 oc |e| £//4 and deduce that 



or oc e M 2 



(5.4.22) 



(see figure |17|) . 



e = All the fixed points have collapsed at the origin around which the motion is stable with 
a positive rotating speed. 



e/a2 > o is the only real fixed point. There are 21 fixed points in the complex plane 



Pk 



dcf 



/ = -e/a 2 < 
\ 9 k = kir/t 



( 



+ 0^(3) 



Pk = i^/e/a2~ cos(k-K / £) 



\ Ik = f i\/ e / a 2 sin(fe7r/£) 



+ 0^(2) (5.4.23) 



for k G {0, . . .,2£ — 1}. 



If we go back from the rotating frame to the fixed one, we add a global rotation of angular 
speed ujq = 2-7rr/(£r). This will make the points p belong to two different periodic orbits of period It. 
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5.5 Summary 

Except the case I = 3 where there are none, in all other cases there are complex periodic orbits in 
the vicinity of the bifurcation. Each of them is connected continuously through the bifurcation with a 
real one. When £ = 4 and a + b < 0, on each side of the bifurcation there is one real and one complex 
satellite orbit around the central point. As the external parameter goes through the bifurcation these 
two orbits interchange, the real one becoming complex and vice versa. In all other cases, when they 
exist, the complex orbits are present on one side of the bifurcation only. 

The bifurcating (real and complex) orbits are typically separated in phase space by a scale e* 
given by Eq. ( [5,4.21| ), This quantity will play a crucial role in the semiclassical considerations of the 
following sections. In the particular case of stable bifurcating orbits, in addition to e* there exists 
another typical scale oc |e|^ 4-1 / 2 which gives the radial width of the stable satellite islands. For the 
extremal (resp. transitional) bifurcation the width of the stable island(s) is given by |e| 3 ^ 4 (resp. |e|). 
In all cases, the classical flow looks roughly the same in both sides of the bifurcation if one looks at it 
with a resolution of order e* or worst. 

6 Normal form for one-dimensional unfolding of generating func- 



So far, we have described the r-periodic dynamics around an equilibrium point by simplifying the 
hamiltonian in its neighborhood. We have revisited the classification of the generic fixed points by 
describing, in each case, the normal form of a hamiltonian whose dynamics, viewed every — 2r, is 
the same as the original one. The one-freedom dynamics thus represents the local Poincare reduction 
around a r-periodic orbit. The associated local Poincare map <&(w;e) is the stroboscopic view of the 
reduced dynamics taken every r. In order to complete our study, in this section we deduce the normal 
forms of the generating function of <&(w;e). In some cases, like in semiclassical theory, the Poincare 
map plays actually a more central role than the reduced hamiltonian itself. 

The starting point for calculating the Poincare map from a time-independent hamiltonian is the 
formula 



tions 




w 



(6.0.1) 



w + t\w, H(w; e)} + J2 —l...{w, H(w; e)j, . . . , H(w; e)\ 



(6.0.2) 



m times. 



47 



where we define the Poisson bracket of two smooth functions on phase space by 

{A,B\ = *(V 1B A)JV„B; (6.0.3) 

drf 8AdB__dB_dA 
dq dp dq dp 

Let idoG^ and construct the backward and forward sequence Vra G Z; w n = f <&( n )(u;o, e)- Some 
relevant functions containing the same information as $(")(• ;e) are 

• The gg-generating function 

F^\w n ,w ) = (p(t)q(t) - H[p(t),q(t);t;e]) dt . (6.0.5) 



w(t) = f ( p fl) is the unique solution of ( 3.0.5| ) such that w(0) = wq and w(nr) = w n . Locally 



depends only on two independent variables among (p n , q n ,Po, Qo)- These are often taken 
to be q n and qo. The values of this function are canonically invariant, they do not depend on 
the choice of the symplectic coordinate chart - possibly time dependent - in which the integral 
in calculated. When H is constant we get 

pnr 

F<ri(w n ,wo) = -nrH(w ;£) + J p{q,H}dt. (6.0.6) 

=H(w n ;e) 

• The mixed generating function 

S {n) (p n ,q ;e) d = p n q n (p n ,q ;e) - F^ipn^o-^) . (6.0.7) 

In all cases but the unstable one, we will not be able to complete the summation ( |6.0.2| ) exactly. 
Nevertheless following our general philosophy we will work at a finite precision. From section [5] we 
have determined for each bifurcation the typical scale in e entering the dynamics. This scale is fixed 
by the order of truncation of the normal forms of the hamiltonian: the more terms are kept, the finer 
details can be described. Accordingly, we will now show that all the dynamical structures up to a 



fixed scale are actually contained in the first terms of the expansion (6.0.2). This allows to compute 
the generating functions for each type of bifurcation. Let us start with the unstable case which is 
special, as we said, in the sense that the sum can be computed explicitly to all orders. The physical 
reason for this is that in a generic unfolding of an unstable point there is no bifurcation and hence no 
local typical scale. 
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6.1 Unstable case 

Since no structural change is induced when e varies, we will leave the e-dependence implicit up to the 
end of this subsection. The one-parameter unfolding is given by ( [4.1.1| ) 

n 

H u (p,q;t) = Xpq + Y, ^k{pq? + M {2n + 1) (6.1.1) 

k=2 



d = K n (pq) 

We recall (see subsection |3.2| ) that the true dynamics is recovered by taking a stroboscopic view 
every \st\, where ? = 1 for the hyperbolic case and ? = ±2 for the inverse hyperbolic case. (In 
the following discussion we take, for the inverse hyperbolic case, the plus sign). Since d(pq)/dt = 
{pq,H u (p,q;t)} £ O m (2n + 1) 

pq is a constant of motion modulo Pj(? (2n + 1) and the solutions of the equations of motion are 
simply 

p(t) = p e- tK - (poqo) + O po , qo (2n) ; (6.1.2a) 
q(t) = q e +tK - (poqo) + O po , qo (2n) . (6.1.2b) 



From ( |6.0.6 ) we get 

Ff?\wi,Wo) = sr[p qo K' n (p q ) - K n (p q )] + O P0j(?0 (2n + 1) . (6.1.3) 

The mixed generating function is obtained from the definition ( |6.0.7l ) after inverting equations ( |6.1.2| ) 
up to first order 

S^\p<:, qo) = PoQo -r[K n (p q ) - p q K' n (p q )] (6.1.4) 

= e' T V?o + ?rh 2 e 2fTA (p,q ) 2 + [?rh 3 + 2( ? rh 2 ) 2 ]e 3?rA (p,q ) 3 + O p ^ go (8). (6.1.5) 

In the inverse hyperbolic case the hamiltonian ( |6.1.1| ) cannot alone give all the information about the 
dynamics at time r. In particular, w(t) given by ( gXp does not contain one transformation whose 
linear part is a symmetry with respect to the origin. If we instead propose as a guess the following 
one-step generating function 

S^{pi,qo) = (Tipiqo + 02{pm) 2 + cr 3 {piq ) 3 + O puqo {8) (6.1.6) 

then the real coefficients a±, 02 and o% can be determined in the following way 
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1- Invert equations ( 3.2. l| ) modulo O pl)9l (7) to find out wi( 



w ) 



, ] , deduce that o\ = — e ; 
o -e / 

3- Iterate the mapping in order to get u>2 modulo O pim {7) and identify the expression with w(2t) 
obtained via equations (|6.1.2|) ; 

4- Eventually find out 02 and 03. 

For the inverse hyperbolic case, we thus obtain 

SPfri,qa) = -e TX piqo ~ rh 2 e 2rX ( Pl q ) 2 + [rh 3 + ^(rh 2 ) 2 ] e 3rA ( Pl q ) 3 + O pi>qo (8) . (6.1.7) 

Note that, in any case, the mixed generating functions have the same functional form as the 
hamiltonian but the relations between their respective coefficients are non-trivial. 

6.2 Extremal case 



From hamiltonian ( 5.2. lj ) the typical scales are p tx e± and q oc \/\e\- We deduce the following Poisson 
brackets 

fotfa} = \p + 0*^(3); (6.2.1) 
fei?x} = -e-aq 2 + O p q ^(3) ; (6.2.2) 
#1} = -^pq + O pf9>v ^|(3) ; (6.2.3) 
l-.-lp, ^},...,^} G O pi9) ^(3). (6.2.4) 

n times, n^3 

From ( |6~a2D we get the Poincare map 



^1 < 



' Pi =Po-T£-TaQo- -ap 2 l -Tap q + ° P0i q 0i% /j7j(3) 5 (6.2.5a) 
ft = qo +Po ~ - \ a( & ~ J^ a Po ~ ^ a Po<lo + ° POi(?0i ^/j7j( 3 ) • (6.2.5b) 



To get the mixed generating function we have two methods 

• The first is to write the general form of S^(pi, on) modulo O mi^) and deduce the Corre- 

ct Po,<?o,\/l £ l 

sponding map. We get S by identifying the latter with ( |6.2.5 ); 
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The second is to calculate S^ 1 ' via the formula ( |3.0.7 ). 



We obtain 



S[ 1] (pi,qo) = ^Pl + ^epi+reqo + Piq + — ap{ + - ap(q + - a Pl q£ + - aq% + (3^ _^(4) , 

(6.2.6) 

and for the gg-generating function ( |6.0.5| ) 



F 



(i) 



1 r 2r t 

Po + re q + rep + - a q% + — aplq + r ap ql + -apl + O 



2 ru " IU ' ■ f " ' 3 ^ u ' 3 *~ u ™ ' ■ • ru ^ u ' 6 ^ u ' po,9o,VI £ I 
Both functions have a functional form quite different from the one of H\. 



(4) . (6.2.7) 



6.3 Transitional case 



From hamiltonian ( 5 . 3 . 1| ) the typical scales are p oc |e| and q oc \/\e\- Remember that this hamiltonian 
describes the true dynamics when we look at the dynamics backwards every 2r. Moreover c; = —2 in 
equation ( 3.1.10| ) and the monodromy matrix for one period is 

(-1 



V 



1 -1 

From the hamiltonian H 2 , we deduce the following Poisson brackets 

\q,H 2 ) 



1 



p+ O yr-r .— (5) ; 



{p, #2} = -2e g - a q 3 - eb q 2 + c q 4 + q> ^-(5) 

{fo^ 2 [,^ 2 [ = -%(2e + 3a( Z 2 )+ 0^^(5) ; 



n times, n^3 

From ( |6.0.2 ) calculated with r 1— > — 2r, we get the Poincare map 



(6.3.1) 



(6.3.2) 
(6.3.3) 
(6.3.4) 
(6.3.5) 



p 2 =p + 4req Q + Iraqi ~ At£ Po + 2rebgg + 2rcq% - 6"rcp <7o + °^|,g , v ^y( 5 ); (6.3.6a) 



(2) 



q 2 = qo - 2p - 4reg - 2raql + -T£ Po - 2rebql - 2rcql + 4-rcp <?o + ^y^i,^, ^( 5 )- 



(6.3.6b) 
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We now compute the one-step mixed generating function. Since it is smooth with respect to e, its 
general form is 



S<P(pi, Qo) = ^VQo + IPiQo + g + 2 Sp i + Kpiq o + 4 01 q o + CzPiQo + ^Piffo + \ L 1o 



1 



1 



+ o AiPi + - A 2 p?g§ + XsPiQo + O 



where (6\, £3, 62, l, X%, A2, A3) £ [R 6 and (77, 7, £, 5, k) are real functions of e 



(6.3.7) 



(5). 



77 = r?o + ; 7 = To + ^71 ; £ = 6 + e£i ; 
(5 = 5o + eSi ; k = ko + e«i ; 

with (?7o,»7i>7o,7l>Co,&>£o>£l,«o,«l) € R 6 and 70 / 



(6.3.? 



Although they are of order 



M>90,\/l£ 



(6), the last three terms in must be taken into ac- 



count. This is because when approximating the Poincare map up to order four (see equations Q4,2.15| ) 
and ( 3.2. 1| ) ) differentiation with respect to p\ makes the order fall by two. Frompo(pii Qo) and gi(pi, go) 



modulo O n — 1 /r-r(5), it is straightforward to invert them to obtain pi(po><7o) an d qi(po,Qo) m( > 
dulo O 



^(5). Condition ( |6~3~TD imposes 

7o = #0 



-1 . 



(6.3.9) 



From this we can calculate P2(po,Qo) and q2(po, Qo) modulo O r, — r /tt(5)- Identifying these two 

V \P0 I 1 Q0 1 "y \£\ 



latter equations with the Poincare map $2 given by ( [4.2.15 ) we get a set of equations with a unique 
solution if and only if 

dh 3 



b = 3— - =0 

as \e=o 



and c = 5h 5 (0) = 0. 



(cf the discussion of subsection 4.2). We finally find 



as well as the Poincare map for one period 



(6.3.10) 



S^iphQo) = -PiQo-reql- - <p\ - - ra q% - rspiqo - - rapiql + °^7f, qo , v ^y( 5 ) ' (6.3.11) 



Pi 



Ql 



-po - 2req - ragg + rep + - reapoql + °y^\, qo , ^wS 5 ^ 



1 3 1 2 
i + po + te q + - ra q$ - - repo + en x q 

+ Ai pl ~ x arpoql + (A 3 - 2(3™) q% + O 



|po[,go, y/W\ 



(5). 



(6.3.12a) 
(6.3.12b) 
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We should have worked to a higher order to determine («i, X%, A3). To complete the study, let us 
give the expression of the two-step mixed generating function 



S?\p2,qo) =P2Q2 - 2tH 2 - [ T p(t)q(t) dt , (6.3.13) 

J 

where p(t) and q(t) are calculated via ( |6.0.2| ). We obtain 

Si 2) (ft, qo) = P2Q0 ~ P\ ~ 2re q 2 - \ r aq A + ^ ^_(5) , (6.3.14) 
as well as the gg-generating function 

FP =si 2 \p 2 ,q )-p 2 q 2 = 2rH 2 (p 2 ,q 2 ;T;e)-2p 2 2 + O VWlq2 ^(5); (6.3.15) 

= -p 2 2 + 2req 2 + -aq 4 + /r-f ^(5). (6.3.16) 

6.4 Stable cases 

We will treat all stable cases for i ^ 3 simultaneously by regrouping hamiltonians ( |5,4.1b| ), ( |5.4.6b| ) 
and ( |4.1.3b| ) in the unified form 

H l {I,0;t;e) = eI + £ a fc (e)/ fc + 6<(e) / £ / 2 cos(^) + 0^^(1 + 1). (6.4.1) 
fee N 



We have seen that both p and q scale as e* given by ( |5.4.21| ). Recall (cf Eq. ( 3.3.13| )) that this normal 



form is obtained looking at the dynamics in a rotating frame of angular speed ujq given by Eq. ( [4.1.21) . 
As in the two previous subsections, let us start computing the Poisson brackets in ( 6.0.2; ) 



\I,H t \ = ih(e) I 1 ' 2 sin(^) + Ojj^e + l); (6.4.2) 
{9, H e \ = e + Y, ka k {e)I k - 1 + \lb t {e) I 1 ' 2 ' 1 cob(£0) + 0^(1-1); (6.4.3) 



ke N 

2^k^£/2 

{...{I,H t },...,H e \ e 0^(1+1); (6.4.4) 

V 

n times, n^2 

{...{9,H £ },...,H £ } G O/j^+1). (6.4.5) 

n times, n^2 

Then, in the rotating frame we obtain I(t) and 9(t) modulo O £j .(£+l) and O ^ £ *(^— 1) respectively. 
Going back to the non-rotating frame by adding uj$t to 9{t) we have 
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I(t) = Jo + £tb e (e)I e o /2 S m[£(0 o + u o t)] + 0^ m (£ + l); (6.4.6a) 

0{t) = 9 + (oo + e)t + ty ka k (e)I*- 1 + ^ t£ b e {e) fj 2 ' 1 cos[£{0 o + u t)] 

fceN 2 ( 6 - 4 - 6b ) 

2 ^ /2 + o^^-i). 



For every integer n, it is then straightforward to express the first orders of 9 n = 6 {jit) and Jo i n 
terms of I n = /(nr)and 9q. The mixed generating function for time n follows immediately 

sf\l n M = I n 6 + n{uj + e)Tl n + nr ^ a k {e) 1^ + nr6 £ (e) I'J 2 cos(£0 o ) 

feeN (6.4.7) 

+ 0^,^ + 1) 

= / n # + nuj Q Tl n + TiTH e {I n ,d ) . (6.4.8) 
We finally calculate the gg-generating function from, for instance, its definition ( |6.0.5| ) 

F i n) = titJ2 (k- l)o fc (e) /o + nr (| " X ) cos (^o) + °v^, + !) • ( 6A9 ) 

fceN V ^ 

2^fc^£/2 

7 Asymptotic expansion of oscillating integrals 

In the previous sections, we have exclusively dealt with classical mechanics. The semiclassical analogue 
of these geometrical studies is highly non-trivial for essentially two reasons. The first is connected to 
the choice of the quantum hamiltonian one must start with. Infinitely many quantum hamiltonians 
differing from each other by /i-dependent terms correspond to the same classical dynamics. The second 
is that the semiclassical analogue of the Poincare reduction (section 0) is very difficult to extend beyond 



the leading order (see [1C, [SJ). Nevertheless, we will assume in the following that one has already 
solved these problems for a given system. 

Let us then begin with a given evolution operator Uj{p, q) depending on one external parameter 7. 
The Hilbert space on which the canonical operators p and q are defined can be finite or infinite. The 
associated classical dynamics corresponds to a one-step iteration Poincare map 

Our aim now is to use Meyer's classification to write down semiclassical expansions of physical 
quantities related to the exact quantum evolution operator U~. For instance, we may compute semi- 

54 



classically the trace of correlation operators of the form 

M 



C{i) = Wa^q)^ (7.0.10) 



def 

i=0 

where M is a strictly positive integer, {a>i} ie {o,...,m} are some smooth functions which do not depend 
on h and {rij}j e {o,...,m} are some integers whose sum is I ^ 0. Without loss of generality we can 
assume £ > by possibly changing the arrow of time. One particular case of such operators are of 
course the trace of the £th power of C/ 7 , which play a central role for calculating the spectral function 



of the quantum system [48 |. 
Formally we have 

QtQa-r) ^ tr[C( 7 )] =E h " [[ VtJktfMi) e^V'^dp'dq , (7.0.11) 

where the F : s are some smooth functions. The integration domain is the whole bidimensional phase 
space and the argument of the exponential is 

N® q . 7 ) = N W (0, 0; 7 ) + p'q - S® (p' , q; 7 ) (7.0.12) 

where is the mixed generating function of the £th iterated Poincare classical map <f>)y = 
Without loss of generality we can take A^(0,0; 7 ) = 0. 

The basic steps leading to Eq. ( |7.0.11|) are the following. First, the fact that #^ (p, q) is the 
classical dynamics associated to Uy(p,q) means that we have 

ffi^ ~ e^V*?* y u W(p',q)h k (7.0.13) 
(P>\q) 

where the u's are some smooth functions except, perhaps, for a zero-measure set. For to be unitary, 
uq has to be the square root of d 2 N^ /dp'dq Appendix A]. 

Second, remark that for every smooth a(p, q) and for every integer n, we have the formal asymptotic 
series 

(p\U- n a(p,q)U%\q) 



7 

~ a 



(p\q) h-*0 



[^ ] (p,q)] + J2aM(p,q)h k (7.0.14) 



for some functions a"'s which are smooth everywhere except for a zero- measure set. This last 
formula can be understood if we write the formal expansion of a(p,q) in powers of (p,q), then use 
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the commutation relation between p and q and eventually make a stationary phase expansion at all 
orders in h after having inserted the appropriate resolution of identity made with the projectors | p ) ( p \ 
or | q)(q |. Except for a zero-measure set of (p, g)'s the stationary points of the oscillating integral are 
non-degenerated. 

Third, write Qi{h]^f) in the following manner 

Qtfa-r) = tT[U^ ai U^a 2 U- ni U^ +n2 a 3 U- n '- n2 ---U^ , (7.0.15) 

and insert alternatively the resolution of identity made with the projectors \p){p\ or \q){q\. We are 
therefore led to a many-dimensional oscillating integral. A stationary phase expansion can be made 
formally to all orders in h in the Morse (quadratic) directions of the critical points. Eq. ( [7.0.1 1| ) is 
thus recovered, where the integral in the quasi-degenerate directions remains to be done 3 . 



If we make a stationary phase approximation of the integrals appearing in ( 7.0.11 ), we see that the 
domains of phase space which contribute are the neighborhoods of the fixed points of which are 
the periodic orbits of whose period divides I. When we want to calculate the contribution to Qe 
of the k repetitions of an mT-periodic orbit {mk = £) we should have specified the normal form of the 



generating functions corresponding to the mapping <J>m (in section |6] it has been done only when 



m 



for extremal and transitional cases). But it is clear, from the series ( |6TQ72D , that chan ging rriT to mkr 



will only modify the coefficients of the polynomials in (p, q, e) but not their form (in the stable case 
we get a simple relation of proportionality: for Fm the r-dependence is given by It, see Eq. ( |6.4.9D ). 
Since we are mainly interested in the ^-dependence of the terms in the semiclassical expansion we 
are allowed to consider only primitive periodic orbits. Therefore we will take for the generating 
functions calculated in the previous section. 

At this stage, we must recall the hierarchy of structures which are present in a generic mixed 
classical dynamics. In a compact domain of phase space there exists a finite number of periodic orbits 
of a given period I. Let us consider one of these orbits. If it is far from its bifurcation point, then 
complex periodic orbits are neglected and the real ones — which are isolated — contribute in the usual 
way |2(| . On the contrary, near its bifurcation one can always identify a typical distance e* which was 
introduced and calculated for each type of bifurcation in ||. Moreover there exists a typical phase-space 
quantum scale A cx \ffi. If A > £*, then a uniform approximation is needed in order to compute the 
semiclassical contribution of the orbits to Qe K0, 49, [5], ^7], 0]. Finally, if A < £*, then one can 



still consider the satellites orbits are isolated but now the complex ones are no longer negligible. This 
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is the typical situation we will deal with in the remaining of this section. More precisely, we will 
evaluate the semiclassical contribution of the real as well as the complex periodic orbits for each type 
of bifurcation described previously. 

After having constructed a suitable symplectic chart, we are therefore led to oscillating integrals 
of the form 

T. / / W(rJ n\A N M*&Ar! 



F(p',q)e^ e{p ' g]£, dp'dq , (7.0.16) 

where F is a smooth generic function. Without loss of generality we can suppose F(0,0) = 1. Ni, de- 
fined by ( [7.0. 12p , is obtained from the normal forms of the mixed generating function sjp of section ||. 
The integration domain is a neighborhood of the fixed island whose central point is at the origin o. It 
is assumed not to overlap with other islands, e must now be considered to be a fixed quantity small 
compared to all other classical quantities. We will control semiclassical errors by varying K only. 

7.1 Extremal case 



From ( fT2l| we have 



-Ni(p,q;e) = T - ep + re q + ^p 2 + ^ap 3 + ^ ap 2 q + T -apq 2 + ^aq 3 + 0,^^(4) ■ (7.1.1) 

Near the fixed points we have F(p,q) = 1 + O p>5 (l) and therefore the uniform approximation 
of ( gjUg) is an Airy function. The two fixed points are isolated if 



h te . (7.1.2) 

They both give an oscillating contribution if they are real (e/a < 0) 

2i = 27dh (e-'fe-^* 1 ^ [1 + 0^(3)] + e 2i ^l 3/2 /^ a l 1/2 ) [l + 0^(3)]) . (7.1.3) 

This is the usual Gutzwiller contributions where we have made explicit the Maslov index for each 
orbit (see figure jig a)). 

If e/a > 0, the integration contour can be deformed to pass through only one complex saddle 
point, giving an exponentially small contribution (see figure |l8| b)) 

T _ 27ri?i _i| _ 2T | e |3/2 /(3ft | |i/2 ) 



v / 2r|ea| 1 /2 



'(l + O ft (2)). (7.1.4) 
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(a) e/a < 



(b) e/a > 



Figure 18: Deformations of the contour through the saddles points for the oscillating integrals X\. 



7.2 Transitional case 
From ( |6.3,14| ) we have 

^2 (P, q)=p 2 + 2re q 2 + ^raq A + g> ^-(5) 

The satellite points can be considered to be isolated from the origin if 

h < T-r-r . 

\a\ 

The contribution of the central fixed point is 



where the Maslov index is given by 



r i def 
V\X\ = < 



if x > 0; 

1 if x < 0. 



If e/a < the satellite points are real and contribute by 



2 / irifi -if ^2We 2 /(ah) 

\ 8r|e| 



(7.2.1) 



(7.2.2) 



(7.2.3) 



(7.2.4) 



(7.2.5) 



5N 



If e/a > we must study the relative orientation of the saddle points. Indeed after having done one 
stationary integration in the p direction, we have 



'Tvih 



+ °° e f(^ 2 W) dq[l + 0vn (2)] 



(7.2.6) 



The saddle points (sq = 0, s±i = iq±) are represented in figure IS. 




forbidden regions 
<3 [> : Directions of steepest descent 



Figure 19: Saddle points for I2 when e/a > 0. 



The hatched regions represent the forbidden regions where 9[-/V(z) — iV(s)] < 0. It is clear from 
the figure that it is impossible to deform the contour ^in order to cross s\ or s_i and to come back 
to the real axis. The complex saddles will therefore contribute to semiclassical expansions only via 
power laws in h when one tries to go beyond the dominant order of the semiclassical expansion at so- 
This power law expansion will contain some information about the presence of s±\. Let us specify the 
corresponding power of h 
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12 = I Li 1 + \ p2d ™ F \^ + \<? d w F \m + ^ «^ 4 ) e^ {2£Tq2+p2) dpdq[l + 0^(3)] 



1 + i m ^^(0,0) + 2^^,9^1(0,0) - ^2 ) + °v^( 3 ) 



(7.2.7) 



It is the third term of the parenthesis which is the explicit leading contribution of the complex satellite 
points. It corresponds to a term h smaller than the leading Gutzwiller contribution of the real one. 
Even if a is not known its sign can be determined by looking at the stability of the central orbit (see 



figure [LlD 



The fact that the complex saddle points are indeed encoded in the asymptotic series of Z2 is well 
known: When F = 1, one can re-sum the complete asymptotic expansion of Z2 and give to it a sense 
in the whole complex /i-plane except over one half line. The discontinuity of I2 between the two sides 
of this cut can be interpreted in terms of an integral whose contour goes actually trough the saddles s\ 



or s_i. See for example |ll|, §2.1 a] or ]2£|, §9.4 a]. 



7.3 Stable cases 

From ( |6.4.7D we have in action-angle variables 



-N e (I,9;e) = 1st I + It ^ a k {e) I k + £Tb e (e) I e/2 cos(^) + O^^t + 1) . (7.3.1) 



fce N 



We are here systematically ignoring problems related to the quantization procedure. In expres- 
sions like ( ffiOlel ) the passage to action-angle variables or to any other symplectic coordinate system 



of jacobian unity can be done. If one starts from the classical problem expressed in action-angle 
coordinates, we may face however some difficulties to quantize it [14| because the phase space has 



a cylindrical topology and therefore cannot be described by a unique symplectic chart. In Feynman 
path integrals or their discretized analogues, this implies to take care of the homotopy of the paths. 
The interfering contribution of each of them is weighted by the character of its homotopy group [p| . 
For the cylinder, it corresponds to a infinite sum over the winding numbers. Semiclassically this will 
select exactly the central orbit whose winding number is r = £tcu/(2tt), the others being at a distance 
much larger than e and a fortiori than h; we are therefore led to ( |7.3.1[ ). For a concrete illustration 
see ||, §2]. 
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The contribution to the central fixed point I = is given by 



2n roc 



F(I,i 



1=0 



1-y £ a k (e)I k + lrbi(e) i"* 



cos 



ke N 



+ 



a k (e) I k + £rb e (e) l l I 2 cos(£9) ) + • 



'ke N 

2<k<£/2 



e~ iiTEl/h dIde . 



(7.3.2) 



The dominant term is —j^- Many powers of h follow among which we want to determine the first 
which depends on bg. If we expand F in Fourier series we obtain 

oo 

F{I,9) = J2 [F k (I)cos(k9)+F k (I) S in(ke)] . 

k=0 



(7.3.3) 



Moreover, we know from the smooth behavior of F in a neighborhood of the origin that we have 
F(I,8) = F ( V2I cos 9,V2I sin^) and thus that, Vk £ IH, F k (I) = A k I k / 2 for some complex con- 
stant A k . In (|7.3.2 ) we get a term of the form 

/2+1 rao 

-^bi(iz) ^ + 0^(21 + 1). (7.3.4) 



-r-Kbi 

n 



h 

Ire 



hx \ - x 

JTeY dX 



But this is not the largest term depending on be since there is a term proportional to ft 



-1 



I 1 * 



2ir poo 
=0 J 1=0 



I* cos 2 Ue)e- iiT£l/H dId9 = -Trb] I - 

IT 



nV- 1 f\ 



-1 • 



(7.3.5) 



When the satellite points are real, then their contribution to semiclassical expansion corresponds 
to the standard Gutzwiller oscillating terms. To leading order one finds 
- For £ = 3 



Z 8 = — [1 + 0^(3)] + e"^ M^/ 8 e- W/(9 ^) [1 + 0^(3)] 



re 



(7.3.6) 



For I = 4 and with the notations of subsection 5.4.2, we can have one or two real families of 



fixed points corresponding to cos (49) = ±1. When they are real their contribution is 
■ 2n[h e A^ a ^+^ a ^)^ 2 e - 2i ^ 2 ^+7rhy h [i + ^ /n (3)] . 



a—b 
a+c 



1/2 



(7.3.7) 



re 
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For i ^ 5 both stable and unstable satellite orbits give the same contribution when eja^ < 0, 



2-irih 



2a? 



e/4 



-i(u[e]+u[±b e ])-K/2 e iere 2 /(4a 2 h) ^ + 0^(3)] 



(7.3.? 



Martin Sieber showed [46| that if one wants to compute the difference between the contribution 
of the stable and unstable satellite orbits as a function of e one needs to work with normal forms up 
to O^q *^(2£ — 4). The expression of this normal form is given by equation (14) of Ref.[46] 

When the satellite orbits are complex, they actually correspond to real negative action I. The 
integration contour in / goes from zero to +oo and therefore cannot cross the saddle points on the 
negative real axis. The contribution of the complex orbits comes entirely from the saddle expansion 
at the origin. It thus contains information about the complex solution through a term of order h 1 



given by Eq. ( 7.3.5 ). 



oc \e/a 2 \ 







SI 



Figure 20: Saddle points for Ii when ejai > and t ^ 4. 



8 Conclusion 



In the first part of this paper we have revisited Meyer's classification of bifurcations in two-freedom 
systems by systematically including and locating the complex periodic orbits which are present on 
each side of the bifurcation (subsection presents a summary of the main features concerning 
this point). Thereafter we have explicitly studied the contribution of complex periodic orbits to 
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semiclassical expansions ( |l.0.l| ) when the classical regime is a mixture of regular and irregular motion. 
Our starting hypothesis was that the ghost orbits which may generically contribute are those which are 
not too deeply located in the complex phase space and therefore can be associated to a bifurcation. We 
have shown that the only case where they do not give a power law contribution in h is the extremal 
bifurcation where the monodromy matrix has a determinant equal to one. Their contribution is 
exponentially small in h (see Q7.1.4 )), whereas the argument of the exponential is proportional to 



the distance from the bifurcation to the power 3/2. Because of this large power and because no real 
orbits come to shadow them, these ghost orbits can be easily observed when one compares exact and 



semiclassical quantities [[3(J, 31]. In the latter reference it has moreover been shown that these orbits 
may also play an important role in the semiclassical description of tunneling effects. In all the other 
types of bifurcations the complex orbits do not create new peaks in the Fourier transform of the 
density of states but rather modify the amplitude of the peak associated to the central real orbit. 

We have systematically ignored in our discussion the consequences a symmetry may have on 
bifurcations. Because discrete as well as continuous symmetries are, from a technical point of view, 
equivalent to the introduction of another continuous parameter in the unfolding in addition to e, their 
presence increases drastically the number of the possible types of bifurcations. One should for example 
consider cases where the linearized motion corresponds to L = (or M = 1), and should therefore 
classify first all the possible cubic hamiltonians. To our knowledge there is no systematic classification 
of the bifurcations of hamiltonians of codimension two or higher. Nevertheless some normal forms for 
certain types of bifurcations can be found in the literature [36| . The simplest cases are of course the 
normal forms given in section ||| with some of the higher order coefficients put to zero. Then we get 
some additional phase portraits like the one obtained in Ref. pj] , table 1 case c]. One can easily go in 
these cases through a semiclassical analysis similar to the one presented here. 

A.M. would like to express his gratitude to Eugene Bogomolny, Stephen Creagh, Martin Sieber and 
Daniel Rouben for many fruitful discussions. He also thanks le Service de Prestations Informatiques du 
Departement de Mathematiques et d'Informatique de l'Bcole Normale Superieure (Paris) for assistance. 
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Notes 

1 In the literature 0, App. 27], one can find parabolic fixed point. 

2 In the literature g App. 27] and §, § 7], one can find elliptic fixed point. 

3 Since we are working with a one-dimensional unfolding there is only one direction in the integration 
domain where a critical point of a smooth function can be non-quadratic and this result is independent 



of the number of arguments of the function [[4^, chap 4, § 5: the splitting lemma]. In ( 7.0.11 ) we 
nevertheless keep two directions of integration for stressing the physical interpretation of the critical 
points in term of the periodic orbits of the dynamics but it can be checked in the following that in all 
cases (extremal, transitional and stable) there is only one direction in which the integration is not a 
gaussian one. 
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Appendix: Summary tables 



monodromy matrix 



e/a < 



e = 



e/a > 



Extremal case 



spM| £=0 ={1,1} 




oc |e| 



3/4 





Hamiltonian normal form: 

1 „2 i _ _ i 1 „ „3 i 1 ), J 



Hi(p,q;t;e) = ±p 2 +eq + ^aq 6 + ± 6g 4 + P;(? ,^=(5) 



Transitional case 




spM| £=0 ={-1,-1} 




a < 




a > 






Hamiltonian normal form: 

1 „2 I „ „2 , 1 „ „4 , Ur- rA 



H 2 (p,q;t;e) = £ p 2 + e q 2 + \ a q A + be g 4 + cq 6 + O^,, ^(7) 



Table 2: Generic bifurcations of extremal and transitional points 
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Stable case: spM| e=0 = {e e } 



e/a < 



e = 



e/a > 






Hamiltonian normal forms: 
H 3 (p,q;t;e) = \ e (p 2 + q 2 ) + ± 6(p 3 - 3 M 2 ) + O p ,,, e (4) 



H 3 (I,9;t;e) = el + ± 6 (2/) 3 / 2 cos(3 0) + 0^(4) 



a) a + 6 > 



/3) a + fo < 









Hamiltonian normal forms: 
fT 4 (p, g; t; e) = \ e {p 2 + q 2 ) + \ a{p A + q A ) + \ bp 2 q 2 + O p>qiVm (S) 



H 4 (I, 9; t; e)=el+ |(3a + 6)/ 2 + |(a - 6)/ 2 cos(4 0) + ^(5) 



Table 3: Generic bifurcations in stable cases for £ G {3,4}. 
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Stable case: spM| e=0 = { e e } 



e/a 2 < 



e = 



e/a 2 > 




i > 5 





Hamiltonian normal forms: 



H t (p, q;t;e) = \e (p 2 + q 2 ) +^a fe (e) (p 2 + ^ 2 ) fc + \h eM (e)\ »[(p + 

+0 Pl9 (^ + l) 



fee N 

2^fc^£/2 



iZ>(J, 0; t; e) = e/ + J> fc (e) / fc + b/(e) I* /2 cos(^) + G^f + 1) 



fee N 

2^fc^£/2 



Table 4: Generic bifurcations in stable cases for 1^5. 
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